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1  INTRODUCTION 


Irregular  wave  kinematics  is  of  fundamental  importance  to  many  coastal  en¬ 
gineering  phenomena.  Sediment  transport  and  scour  are  two  processes  highly 
dependent  on  a  good  understanding  of  fluid  particle  velocities.  Incipient  motion 
criteria  are  usually  a  function  of  fluid  velocity.  As  well,  the  phenomena  of  sed¬ 
iment  transport  under  wave  action  alone  is  believed  to  be  partly  a  result  of  the 
nonsymmetrical  shape  of  the  wave  profile  and  associated  nonlinear  velocities.  In 
this  thesis  the  nonlinear  effects  of  the  nonsymmetrical  wave  profile  on  the  water 
particle  velocities  are  investigated. 

One  objective  of  this  research  is  to  develop  empirical  equations  for  determining 
the  parameters  of  probability  distribution  functions  for  particle  velocities.  The 
probability  distribution  functions  for  both  peak  velocities  and  for  the  entire  velocity 
time  series  are  investigated.  Also,  all  analyses  are  completed  for  both  the  horizontal 
and  vertical  components  of  velocity. 

Another  objective  is  to  investigate  a  purely  theoretical  expression  for  the  root- 
mean-squared  velocity,  (urm„  and  wrma) ,  of  the  entire  velocity  time  series.  The 
expression  is  based  on  the  incident  wave  height  spectrum,  reflection  coefficients, 
and  linear  wave  theory.  The  range  of  validity  and  the  effect  of  nonlinearities  are 
studied.  The  expressions  are  simply  the  standard  deviation  of  the  velocity  and  are 
useful  parameters  in  studies  of  turbulence  and  sediment  transport. 

It  is  intended  that  this  research  provide  a  firm  foundation  for  further  studies 
related  to  sediment  transport  and  seafloor  scour.  For  this  reason,  the  experiments 
have  been  performed  in  the  lower  half  of  the  water  column  in  a  tank  with  a 
horizontal  bed.  Thus,  the  results  reported  here  do  not  necessarily  apply  to  wave 
kinematics  near  the  free  surface. 


1.1  Present  Status  of  the  Question 


Previous  investigations  in  wave  induced  fluid  velocities  show  generally  “good” 
agreement  with  linear  theory.  Dean  (1990)  reviews  the  state  of  the  art  in  water 
wave  kinematics.  For  regular  (monochromatic)  waves  it  appears  that  current  theo¬ 
ries  are  adequate  to  within  the  limit  of  errors  inherent  in  the  laboratory  equipment. 
This  limit  is  suggested  by  Dean  to  be  approximately  5%.  For  irregular  waves  the 

Citations  in  the  thesis  will  be  consistent  with  those  of  the  Journal  of  Waterway,  Port,  Coastal 
and  Ocean  Engineering,  ASCE. 
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existing  literature  is  somewhat  more  vague.  Dean  ( 1990)  reports  that  . .  available 
wave  theories  seem  to  represent  adequately  the  measurements”  for  the  study  of 
kinematics  below  the  trough  elevation. 

The  previous  investigations  described  below  summarize  the  literature  most 
relevant  to  this  thesis.  Surprisingly,  much  of  the  work  in  this  thesis  has  not  been 
done  before  despite  being  of  a  fundamental  nature. 

This  thesis  further  develops  the  ideas  in  the  paper  Estimating  Wave-Induced 
Bottom  Velocities  at  a  Vertical  Wall  (Hughes,  1992).  In  this  paper,  the  general 
theoretical  expression  for  the  root-mean-squared  velocity  of  the  entire  time  series 
was  simplified  to  the  case  of  irregular  waves  in  front  of  a  vertical  wall.  The  expres¬ 
sion  was  validated  using  laboratorv  measurements.  However,  no  range  of  validity 
or  empirical  augmentation  for  nonlmearities  was  determined.  Also,  the  paper  only 
investigated  horizontal  velocities  and  assumed  a  perfectly  reflecting  vertical  wall. 

One  paper  closely  related  to  this  investigation  was  that  by  Soulsby  (1987) 
where  a  method  was  presented  for  calculating  the  root-mean-squared  bottom  or¬ 
bital  velocity  given  the  wave  height,  period  and  water  depth  and  any  one  of  the 
following  spectral  forms;  JONSWAP,  Pierson-Moskowitz,  Bretschneider,  ISSC  or 
ITTC.  The  method  of  calculating  urrna  is  identical  to  that  by  Hughes,  except  it 
does  not  allow  for  wave  reflections.  This  paper  is  strictly  theoretical  and  no  exper¬ 
iments  were  performed.  Thus,  the  limits  of  nonlinearities  were  not  studied,  and 
no  experimental  data  were  presented  to  illustrate  the  validity  of  the  theoretical 
expression.  The  paper  only  considered  horizontal  near-bottom  velocities. 


1.1.1  Velocity  Probability  Distribution 


Longuet-Higgins  (1963)  derived  the  probability  density  function  of  the  water 
surface  elevation  using  the  Gram-Charlier  series  to  expand  the  Gaussian  probabil¬ 
ity  density  function.  This  thesis  applies  this  expanded  distribution  to  fluid  particle 
velocities  instead  of  wave  elevations.  Ochi  (1982)  presents  a  review  of  recent  ad¬ 
vances  in  the  stochastic  analysis  of  random  seas.  He  notes  that  the  first  term 
of  the  extended  distribution,  henceforth  called  the  Gram-Charlier  distribution, 
represents  the  Gaussian  waves  and  the  higher  order  terms  are  attributed  to  the 
nonlinear.  non-Gaussian  wave  characteristics.  The  author  notes  that  there  does 
not  seem  to  be  any  strict  theoretical  proof  that  the  higher  order  terms  in  the  series 
are  necessarily  a  result  of  wave  nonlinearities.  It  seems  that  this  explanation  is  de¬ 
rived  through  a  process  of  elimination.  That  is,  it  can  be  proved  mathematically, 
(Longuet-Higgins  1963).  that  purely  linear  waves  will  have  a  Gaussian  distribu¬ 
tion  Therefore,  any  deviation  from  a  Gaussian  distribution  must  be  attributed  to 
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wave  nonlinearities. 


Ochi  (1982)  discusses  a  series  of  experiments  by  Honda  and  Mitsuyasu  (1976) 
comparing  the  Gram-Charlier  distribution  with  measurements  of  wave  elevations, 
particle  velocities  and  particle  accelerations  for  the  entire  time  series.  (Unfortu¬ 
nately,  the  original  paper  is  a  report  from  a  Japanese  university  and  not  available 
to  the  author.)  Reportedly,  very  good  agreement  is  found  for  the  case  of  wave 
elevations  but  not  for  wave  accelerations.  Given  the  close  relationship  between 
acceleration  and  velocity  it  is  reasonable  to  suppose  that  poor  agreement  was  also 
found  for  the  case  of  wave  velocities.  Poor  agreement  might  be  attributable  to  the 
quality  of  measurement  techniques  in  1975. 

Ochi  (1982)  also  covers  two  purely  theoretical  works  concerning  the  probability 
distribution  of  non-Gaussian  wave  elevations.  Tayfun  (1980)  uses  second  order 
Stokes  theory  and  the  assumption  of  a  narrow  banded  wave  spectrum  to  derive  a 
function  that  requires  numerical  integration  for  its  solution.  The  distribution  is 
compared  to  simulated  data.  St.  Denis  (1975)  uses  third  order  Stokes  theory  and 
a  joint  probability  distribution  function  of  wave  heights  and  wave  periods.  The 
wave  heights  and  periods  are  selected  at  random  from  these  two  distributions.  Both 
Tayfun’s  and  St.  Denis’s  methods  may  be  applicable  to  wave  particle  velocities  but 
such  an  investigation  is  beyond  the  scope  of  this  thesis. 

Anastasiou,  Tickell  and  Chaplin  (1982)  reported  that  the  statistical  distribu¬ 
tion  of  both  horizontal  and  vertical  particle  velocities  follow  a  Gaussian  form  for 
points  measured  more  than  half  a  significant  wave  height  below  the  mean  water 
level.  The  measurements  of  velocities  were  made  with  a  laser  Doppler  velocimeter 
(LDV)  under  irregular  waves  in  a  wave  tank.  However,  their  study  focused  on 
velocities  in  the  top  half  of  the  water  column  and  in  the  wave  crest. 


1.1.2  Peak  Velocities 


The  literature  concerning  the  probability  distribution  of  time  series  peaks  (as 
opposed  to  the  entire  time  series)  is  more  extensive.  The  Beta-Rayleigh  distribu¬ 
tion  applied  in  this  thesis  to  velocity  time  series  peaks  was  originally  developed 
for  shallow  water  wave  heights  (Hughes  and  Borgman  1987  and  Borgman  1986). 
Hughes  and  Borgman  reported  a  “reasonable  fit”  to  measured  wave  height  field 
data. 

Ochi  (1982)  and  Hughes  and  Borgman  (1987)  discuss  several  proposed  wave 
height  distributions.  Studies  have  been  done  that  use  the  Gluhovski  distribution, 
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the  Ibrageemov  distribution  and  the  Weibull  distribution.  A  number  of  studies 
have  tried  modifications  to  the  Rayleigh  distribution  and  the  Gamma  distribution. 
However,  all  the  theoretical  studies  look  exclusively  at  the  distribution  of  wave 
heights  rather  than  peak  velocities.  The  author  is  unaware  of  any  attempts  at 
fitting  a  probability  distribution  to  peak  particle  velocities. 

Vis  (1980)  compares  theoretical  and  measured  orbital  velocities  for  irregular 
waves  in  a  wave  tank.  Measurements  were  made  with  a  LDV  and  compared  with 
predicted  velocities  based  on  linear  wave  theory  and  the  synoptically  measured 
water  level  elevations.  Both  horizontal  and  vertical  components  of  velocity  were 
analysed.  An  important  observation  was  that  measured  horizontal  velocities  un¬ 
der  the  trough  appeared  greater  in  magnitude  than  horizontal  velocities  under 
the  crest.  Linear  theory  predicts  symmetrical  velocities  under  the  crest  (particle 
motion  in  the  positive  direction)  and  under  the  trough  (negative  direction).  The 
measured  vertical  velocities  appeared  to  be  more  symmetrical.  Also  noted  was  the 
somewhat  surprising  result  that  linear  wave  theory  more  accurately  predicted  ve¬ 
locities  near  the  free  surface  (just  beneath  the  deepest  trough)  than  at  mid-depth. 
In  addition,  linear  theory  more  accurately  predicted  bottom  velocities  than  those 
at  mid-depth,  as  expected.  Vis  also  notes  that  peak  velocities  appeared  to  be 
nearly  Rayleigh  distributed. 

Daemrich,  et.  al.  (1980,  1982)  also  compared  theoretical  and  measured  orbital 
velocities  for  irregular  waves  in  a  wave  tank.  However,  these  measurements  were 
made  with  inductive  type  and  impeller  type  probes  rather  than  with  a  LDV.  Both 
horizontal  and  vertical  peak  velocities  were  analysed.  This  study  observed  that 
measured  peak  velocities  under  the  trough  were  higher  than  measured  peak  veloc¬ 
ities  under  the  crest,  as  was  observed  by  Vis  (1980).  Neither  Vis  nor  Daemrich, 
et.  al.  quantified  this  discrepency. 


1.1,3  Other  Investigations 


Lee,  Skjelbreia  and  Raichlen  (1982)  compared  theoretical  and  measured  water 
particle  velocities  for  solitary  waves  in  a  wave  tank.  Measurements  were  made 
with  a  LDV  for  both  the  horizontal  and  vertical  components  of  velocity.  Existing 
theories  of  varying  order  were  found  to  agree  well  with  the  measurements. 

Battjes  and  Van  Heteren  (1984)  compared  field  measurements  of  fluid  particle 
velocities  with  linear  wave  theory.  The  velocity  meter  measured  all  three  compo¬ 
nents  of  velocity  using  acoustic  pulses.  Water  surface  elevations  were  measured 
with  a  resistance  type  step  gage.  Possible  errors  in  the  velocity  measurements 
were  reported  not  to  exceed  ±5%.  Measured  vertical  velocities  agreed  very  well 
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with  linear  theory  predictions.  Measured  horizontal  velocities  were  systematically 
lower  than  predicted  velocities  by  an  estimated  5%,  which  is  within  the  range  of 
possible  measurement  error.  Measurements  of  water  particle  velocities  and  water 
surface  elevations  were  not  synoptic. 

Bullock  and  Short  (1985)  compared  theoretical  and  measured  water  particle 
velocities  for  regular  (monochromatic)  waves  in  a  wave  tank.  Both  the  horizontal 
and  vertical  components  of  water  particle  velocity  were  measured  with  a  LDV. 
(The  two  components  could  not  be  measured  simultaneously  with  this  particular 
LDV.)  Comparison  with  Stokes’  theories  of  various  orders  revealed  that  water 
particle  velocities  were  not  always  well  predicted.  Errors  of  over  10%  are  possible  in 
the  magnitude  of  the  first  harmonic  of  velocity  and  greater  errors  may  be  expected 
for  higher  harmonics.  A  significant  mean  horizontal  velocity  was  measured  with 
a  magnitude  of  up  to  20%  of  the  first  harmonic  velocity.  The  investigation  also 
noted  the  importance  of  laboratory  effects  such  as  reflection  from  the  spending 
beach  and  the  existence  of  standing  waves  in  the  wave  tank. 

Swan  (1990)  measured  mean  horizontal  velocities  (also  known  as  drift  velocity 
or  convection)  with  a  LDV  in  a  wave  tank.  The  development  of  a  drift  velocity 
over  time  was  measured.  In  the  early  stages  the  drift  velocity  was  very  similar  to 
Stokes’  solution  for  irrotational  Lagrangian  mass  transport.  As  time  progressed  the 
magnitude  of  the  backflow  steadily  increased  until  the  velocities  were  three  to  four 
times  larger  than  the  initial  mean  velocity.  The  time  scale  for  this  development 
was  on  the  order  of  one  to  ten  hours  and  is  therefore  not  of  immediate  relevance  to 
this  thesis.  However,  the  investigation  did  emphasize  the  importance  of  laboratory 
effects. 

Hansen  (1990)  analysed  existing  experimental  data  on  regular  waves  in  the  surf 
zone  (wave  breaking  conditions).  He  noted  that  measured  horizontal  water  particle 
velocities  are  virtually  independent  of  elevation  in  the  surf  zone.  Also,  linear  theory 
greatly  overestimates  the  magnitude  of  particle  velocity  for  this  condition. 

Kim,  Randall,  Krafft,  and  Boo  (1990)  measured  the  kinematics  of  large  tran¬ 
sient  waves  with  a  LDV.  The  measured  velocities  under  the  crest  were  compared 
with  theoretical  predictions  based  on  an  “equivalent”  Stokes  wave.  It  was  con¬ 
cluded  that  measured  transient  wave  velocities  were  much  more  severe  than  pre¬ 
dicted  velocities  above  the  still  water  level  and  somewhat  less  severe  below  the 
still  water  level. 

Zhang,  Randall,  and  Spell  (1991)  conducted  a  more  general  study  of  wave  crest 
kinematics.  Velocities  were  measured  with  a  LDV  and  compared  with  theoreti¬ 
cal  values  based  on  the  following  three  common  approximate  methods;  Wheeler 
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stretching,  linear  extrapolation,  and  vertical  extrapolation.  It  was  concluded  that 
the  prediction  accuracy  of  each  method  depends  on  the  bandwidth  of  the  wave 
spectrum.  For  narrow  bandwidths,  linear  extrapolation  is  a  better  predictor.  For 
broad  bandwidths,  Wheeler  stretching  is  a  better  predictor. 
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2  EXPERIMENTAL  PROCEDURES 


Laboratory  experiments  were  performed  in  a  flat  bottom  two-dimensional 
wave  tank  at  the  Coastal  Engineering  Research  Center  in  Vicksburg,  Mississippi. 
The  dimensions  of  the  tank  are  46  cm  wide  by  approximately  46  m  long,  (Figures 
la  and  3).  The  height  of  the  walls  is  90  cm.  A  gently  sloping  beach  covered 
with  wave  absorbing  material  was  at  one  end  of  the  tank  opposite  the  wave-maker 
to  minimize  wave  reflections.  Three  layers  of  horsehair  mats  were  used  on  the 
spending  beach  with  a  total  thickness  of  approximately  7  cm.  The  slope  of  the 
beach  was  1:30. 

An  array  of  three  surface  piercing  capacitance  type  wave  gauges  were  placed  in 
the  middle  of  the  wave  tank,  (Figures  2a  and  3).  Each  gauge  measured  the  water 
surface  elevation  as  a  function  of  time.  Incident  and  reflected  wave  height  spectra 
were  calculated  from  these  measurements  using  the  method  of  Goda  and  Suzuki 
(1976).  The  wave  gauges  were  calibrated  before  each  day’s  measurements.  The 
estimated  error  for  the  gauge  calibration  is  ±0.1  cm  (Briggs  1991).  Other  sources 
of  error  in  the  measurement  of  water  elevations  are  due  to  the  internal  electronics 
of  the  VAX  computer  and  are  assumed  negligible. 

Water  particle  velocities  were  obtained  using  a  laser  Doppler  velocimeter  (LDV), 
which  permits  non-intrusive  measurements  of  fluid  velocities.  The  LDV  is  pow¬ 
ered  by  an  Innova  70-2  continuous  Argon-Ion  laser  manufactured  by  Coherent, 
Inc.,  Palo  Alto,  California.  The  wavelength  of  the  emitted  light  is  488.0  nm  and 
514.5  nm,  one  for  each  component  of  velocity  being  measured.  The  laser  permits 
high  sampling  rates  for  velocities  in  both  the  horizontal  and  vertical  directions. 
The  maximum  power  level  is  2  W.  The  output  beam  diameter  is  1.5  mm  for  514.5 
nm  wavelength.  The  laser  is  mounted  on  an  optical  bench  and  can  be  seen  in 
the  lower  left  corner  of  Figure  lb.  Also  visible  in  Figure  lb  is  the  microcomputer 
used  to  record  and  analyze  the  water  particle  velocity  time  series.  The  actual  laser 
beams  are  visible  in  the  lower  T.ght  corners  of  Figures  2a  and  2b  as  the  beams 
emerge  from  the  fiber  optic  probe.  Velocities  were  measured  at  the  laser  beam 
crossing  point,  which  for  these  experiments  was  located  approximately  7  cm  from 
the  inside  face  of  the  glass  side  wall.  Seven  centimeters  from  the  face  of  wall  is  well 
outside  the  boundary  layer  since  for  steady  flow  on  a  smooth  surface  the  boundary 
layer  thickness  is  on  the  order  of  1  mm  (Munson,  Young,  and  Okiishi  1984). 

The  physics  of  laser  doppler  velocimetry  has  been  well  documented  in  other 
sources  (Goldstein  1983).  In  simple  terms,  fluid  particle  velocities  are  measured 
when  a  small  particle  suspended  in  a  fluid  passes  through  the  fixed  volume  defined 
by  the  intersecting  laser  beams.  The  LDV  used  in  this  investigation  has  four  laser 
beams  emitted  by  the  probe,  two  in  the  horizontal  direction  and  two  in  the  vertical 
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(b)  Laser  and  Micro-computer 


figure  1:  Wave  Tank  Photographs,  (a)  Oblique  View;  (b)  Laser  and  Micro¬ 
computer 


(b)  Laser  Probe 


Figure  2:  Equipment  Close-up  Photographs  (a)  Capacitance  Type  Wave 
and  Laser  Probe  (b)  Laser  Probe 
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direction.  Each  laser  beam  has  a  precise  frequency  and  all  the  light  waves  are  in 
phase.  The  intersection  of  the  two  beams,  therefore,  sets  up  an  interference  pattern 
of  alternating  light  and  dark  bands  called  a  “fringe”  pattern.  As  a  particle  passes 
through  this  fringe  pattern  light  is  scattered  when  the  particle  is  in  a  “light”  band 
and  is  not  scattered  when  it  is  in  a  dark  band.  The  spacing  of  the  fringe  pattern 
is  known  from  the  physics  of  light.  By  measuring  the  frequency  of  the  flashes  of 
scattered  light  it  is  in  principle  possible  to  measure  the  velocity  of  the  particle  as 
it  passes  through  the  fringe  pattern.  By  adjusting  the  relative  phase  of  the  two 
intersecting  laser  beams  it  is  possible  to  get  the  fringe  pattern  itself  moving  at  a 
fixed  velocity  thereby  permitting  differentiation  between  velocities  in  the  positive 
and  negative  directions.  The  LDV  used  in  this  investigation  detects  the  flashes 
of  scattered  light  at  the  same  probe  that  emits  the  laser  beams.  This  is  called 
measuring  in  the  backscattering  mode.  The  probe  is  connected  to  the  laser  and 
signal  processing  electronics  by  a  fiber  optic  cable. 

The  doppler  frequency  range  used  was  33  to  333  khz,  which  provides  an  accu¬ 
racy  in  velocity  measurements  of  ±0.015  m/s.  In  addition,  the  internal  electronics 
have  a  constant  offset  that  must  be  corrected.  The  velocity  offsets  are  a  constant 
source  of  error  and  are  measured  by  placing  the  laser  measuring  point  on  the  glass 
sidewalls  of  the  wave  tank.  The  velocity  of  the  glass  is  certainly  zero,  and  there¬ 
fore,  the  velocity  that  the  LDV  measures  can  be  applied  as  an  offset  to  all  future 
measurements.  Typical  offset  values  are  shown  in  Figure  4a  (Figure  4a  is  discussed 
later  in  this  thesis).  Offset  values  were  determined  before  each  days  measurements. 

Titanium  dioxide  particles  from  white  latex  paint  were  used  to  “seed”  the  wave 
tank  in  order  to  provide  particles  for  scattering  the  light.  The  amount  of  seeding 
added  to  the  tank  was  dictated  by  experience  and  a  feel  for  the  quality  of  the 
signal  return.  Too  much  seeding  clouded  the  water  excessively  and  attenuated  the 
laser  beam.  Too  little  seeding  resulted  in  “dropouts”  whereby  no  particle  passed 
through  the  laser  intersection  point  resulting  in  no  signal.  About  one  teaspoon  of 
white  latex  paint  (without  the  binding  agent)  added  to  fresh  tank  water  resulted 
in  a  good  signal  (approximately  1  teaspoon  per  5  m3).  Chlorine  was  also  added 
to  the  tank  to  help  control  algae  growth. 

The  measurements  of  velocities  with  the  LDV  and  of  wave  heights  with  the 
capacitance  type  gauges  were  not  synoptic  That  is.  a  water  particle  velocity  mea¬ 
surement  at  a  certain  instant  in  time  cannot  be  matched  with  a  water  surface 
elevation  measurement  for  the  same  instant  in  time.  LDV  velocities  were  recorded 
and  analyzed  on  a  CompuAdd  386/20MHz  microcomputer  while  water  surface  ele¬ 
vations  were  measured  and  analyzed  on  a  completely  separate  VAX  minicomputer. 
Analysis  of  the  calculated  statistics  for  each  time  series  was  done  mainly  with  the 
aid  of  a  personal  computer  based  spreadsheet.  That  is,  statistics  calculated  from 
the  LDV  measurements  and  from  the  separate  wave  gauge  measurements  were 
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Figure  4:  Computer  Output  -  Water  Particle  Velocity  Time  Series  Analysis 
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Figure  4:  Continued 


Table  1:  Experimental  Plan 


X- coordinate 

-10-5  m 

Z-coordinate 

mid-depth,  near  bottom 

Water  Depth 

50,  30,  20,  15  cm 

Peak  Spectral  Period 

1.0,  1.5,  2.0,  2.5,  3.0  seconds 

Significant  Wave  Height 

4.0,  9.0  cm 

entered  by  hand  onto  a  spreadsheet.  The  spreadsheet  was  mainly  used  for  calcu¬ 
lating  nondimensional  parameters  and  tabulating  data  points.  PC-based  graphics 
software  was  used  for  graphing  the  results. 

Three  wave  gauges  were  used  to  determine  the  incident  and  reflected  wave 
height  spectrums  using  the  method  of  Goda  and  Suzuki  (1976).  The  reflection 
coefficient  is  calculated  by  dividing  the  total  energy  in  the  reflected  spectrum 
by  the  total  energy  in  the  incident  spectrum.  Figure  5  graphs  the  incident  and 
reflected  spectrums.  The  discrete  spectral  components  graphed  in  Figure  5  are 
written  to  a  separate  data  file  for  use  in  calculating  uTm,  and  wrms. 


2.1  Experimental  Plan 


At  the  core  of  the  experimental  plan  are  80  time  series  distinguished  by 
combining  every  way  possible  the  values  of  the  five  variables  listed  in  Table  1. 
Table  2  lists  all  the  time  series  which  were  run  for  irregular  waves. 

Wave  parameters  were  selected  to  cover  a  complete  range  of  wave  conditions 
consistent  with  the  limitations  of  the  wave  tank.  Of  particular  interest  were  non¬ 
linear  wave  conditions,  and  the  experimental  program  was  structured  so  wave 
nonlinearities  more  or  less  gradually  increased  as  the  testing  progressed.  This  the¬ 
sis  is  intended  to  provide  a  basis  for  further  studies  related  to  sediment  transport 
and  scour.  Therefore  the  measurement  points  are  in  all  cases  at  mid-depth  or  near 
bottom.  Table  2  also  lists  several  time  series  that  have  identical  wave  conditions  to 
time  series  run  previously.  These  series  were  run  to  show  repeatability  and  spatial 
homogeneity  in  the  velocity  field 

Not  listed  in  Table  2  are  the  eighteen  regular  wave  time  series  that  were  run 
before  the  irregular  wave  conditions.  These  regular  waves  were  run  in  order  to 
optimize  the  design  of  the  wave  absorbing  spending  beach. 


15 


Spectral  Density  (ft*/Hz) 


Figure  5:  Measured  Incident  and  Reflected  Wave  Spectra 
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Table  2:  Summary  -  Individual  Time  Series 


Measuring  Point 


X-coord.  -  Z-Coord. 
(m)  (cm) 


Planned  Conditions 


7.00 
7.00 
0.50 
0.50 
0.50 
0.50 
7.00 
7.00 
0.50 
0.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-7.00 
-8.75 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 
-10.50 


-10. 
-10. 
-10. 
-10. 
-20. 
-20. 
-20. 
-20. 
-25. 
-25. 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-25.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-47.00 
-15.00 


20.00 

20.00 

20.00 

20.00 

40.00 

40.00 

40.00 

40.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 


Notes 


4.00 
9.00 
4.00 
9.00 
4.00 
9.00 
4.00 
9.00 
4.00 
9.00 
4.00 
9.00 
4.00 
9.00 
9.00 
9.00 
4.00 
4.00 
4.00 
4.00 
4.00 
9.00 
4.00 

4.00  repeat 

9.00 

4.00 

9.00 

4.00 

9.00 

4.00 

9.00 

4.00 

9.00 

4.00 

9.00 

4.00 
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Table  2:  Continued 


Measuring  Point 

Planned  Conditions 

Run 

X- coord. 

Z-Coord. 

im 

mm 

No. 

(m) 

(cm) 

1331 

ISSl 

Notes 

30 

-10.50 

-15.00 

30.00 

1.00 

9.00 

■■ 

31 

-10.50 

-15.00 

30.00 

1.50 

4.00 

32 

-10.50 

-15.00 

30.00 

1.50 

9.00 

33 

-10.50 

-15.00 

30.00 

2.00 

4.00 

mm 

34 

-10.50 

-15.00 

30.00 

2.00 

9.00 

mm 

340 

-7.00 

-15.00 

30.00 

2.00 

9.00 

Em 

341 

-8.75 

-15.00 

30.00 

2.00 

9.00 

35 

-10.50 

-15.00 

30.00 

2.50 

4.00 

36 

-10.50 

-15.00 

30.00 

2.50 

9.00 

37 

-10.50 

-15.00 

30.00 

3.00 

4.00 

38 

-10.50 

-15.00 

30.00 

3.00 

9.00 

39 

-10.50 

-27.00 

30.00 

1.00 

4.00 

40 

-10.50 

-27.00 

30.00 

1.00 

9.00 

41 

-10.50 

-27.00 

30.00 

1.50 

4.00 

42 

-10.50 

-27.00 

30.00 

1.50 

9.00 

43 

-10.50 

-27.00 

30.00 

2.00 

4.00 

44 

-10.50 

-27.00 

30.00 

2.00 

9.00 

45 

-10.50 

-27.00 

30.00 

2.50 

4.00 

46 

-10.50 

-27.00 

30.00 

2.50 

9.00 

47 

-10.50 

-27.00 

30.00 

3.00 

4.00 

48 

-10.50 

-27.00 

30.00 

3.00 

9.00 

49 

-10.50 

-10.00 

20.00 

1.00 

4.00 

50 

-10.50 

-10.00 

20.00 

1.00 

9.00 

51 

-10.50 

-10.00 

20.00 

1.50 

4.00 

52 

-10.50 

-10.00 

20.00 

1.50 

9.00 

53 

-10.50 

-10.00 

20.00 

2.00 

4.00 

54 

-10.50 

-10.00 

20.00 

2.00 

9.00 

55 

-10.50 

-10.00 

20.00 

2.50 

4.00 

56 

-10.50 

-10.00 

20.00 

2.50 

9.00 

57 

-10.50 

-10.00 

20.00 

3.00 

4.00 

58 

-10.50 

-10.00 

20.00 

3.00 

9.00 

580 

-10.50 

-10.00 

20.00 

3.00 

9.00 

repeat 

581 

-10.50 

-10.00 

20.00 

3.00 

9.00 

repeat 

582 

-8.75 

-10.00 

20.00 

3.00 

9.00 

repeat 

583 

-7.00 

-10.00 

20.00 

3.00 

9.00 

repeat 

59 

-10.50 

-17.00 

20.00 

1.00 

4.00 
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Table  2:  Continued 


Measuring  Point 

X-coord. 

Z-coord. 

(m) 

(cm) 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-17.00 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-8.75 

-7.50 

-7.00 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-7.50 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

-10.50 

-12.00 

Planned  Conditions 


20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

20.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 

15.00 
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Each  irregular  wave  time  series  was  run  in  the  wave  tank  for  approximately  5 
min.  Collection  of  velocity  data  with  the  LDV  began  20  s  after  the  start  of  wave 
generation  to  allow  time  for  the  first  waves  to  travel  down  the  tank  and  be  reflected 
back  past  the  wave  gauge  array.  Velocity  data  were  collected  at  a  30  Hz  rate  for 
each  velocity  component  for  273  s.  Collection  of  water  elevation  data  with  the 
wave  gauges  was  at  a  20  H  rate  for  250  s.  The  5  min.  length  of  each  time  series 
was  selected  based  on  experience.  A  long  time  series  is  desired  for  ensuring  stable 
statistics.  However,  a  short  time  series  minimizes  the  buildup  of  reflected  waves. 


2.2  Computer  Output 


Figures  4  and  5  are  hardcopy  outputs  for  one  time  series.  Figure  4  is  for 
velocity  data  measured  by  the  LDV.  Figure  5  is  from  the  VAX  and  is  the  result 
of  the  analysis  of  water  surface  elevations  (tabulated  statistics  from  the  VAX 
computer  are  not  given).  The  calculated  statistics  from  each  of  these  two  outputs 
were  then  loaded  on  a  spreadsheet  for  further  study  of  trends  between  different 
time  series.  (In  Figure  4a,  a  value  of  -0.0404  cm/s  for  the  “U-direction  velocity 
minimum”  refers  to  the  smallest  absolute  value  velocity  measured  in  the  time 
series.) 

Figure  6  shows  a  portion  of  one  measured  velocity  time  series  produced  by  the 
LDV  and  subjected  to  various  filtering  options.  Similar  graphs  were  produced  for 
every  time  series  during  the  experiments  and  examined  visually  on  the  computer 
monitor.  The  plots  provide  a  useful  check  on  the  quality  of  the  velocity  time  series. 
Evident  in  the  raw  data  (Figure  6a)  are  high  frequency  oscillations  caused  by 
turbulence  and  imperfections  in  the  electronic  data  acqT  isition  systems.  The  effect 
of  filtering  the  data  is  clearly  seen  in  Figures  6t  to  6d.  Filtering  is  discusssed  later 
in  this  thesis.  Figure  6a  is  not  typical  of  all  the  time  series.  Approximately  half 
the  runs  looked  like  Figure  6b  pnorto  any  filtering.  Based  on  a  visual  examination 
of  portions  of  all  the  velocity  time  series,  the  author  is  confident  that  the  velocity 
measurements  were  of  a  very  high  quality. 


2.3  Wave  Signal  Generation 


Waves  were  created  in  the  tank  by  a  vertical  board  moved  through  the  water 
by  a  computer  controlled  hydraulic  ram.  The  exact  horizontal  position  of  the 
board  was  updated  every  1/60  of  a  second.  The  board  signal  was  generated  by  a 
VAX  Fortran  program  called  TSGMN3P1.FOR.  The  user  enters  into  the  program 
the  following  parameters: 
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6.0E-Z 


-6 .0E-2-I - ( - 1 - 1 - 1 - ( - 1 - (— 

138  139  148  Ml  MZ  M3  M4  M5 

TINE  -  seconds 


(a)  Raw  Data 


Figure  6:  Velocity  Time  Series  -  Run  No.  19.  (a)  Raw  Data;  (b)  Filtered  Data 
[0.20-5.00j  s;  (c)  Filtered  Data  [0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 


21 


VELOCITY  -  neters/second  VELOCITY  -  weters/'second 


(c)  Filtered  Data  [0.20-10.00]  s 


(d)  Filtered  Data  [0.10-5.00]  s 
Figure  6:  Continued 
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•  Water  Depth 

•  Peak  Spectral  Period 

•  Significant  Wave  Height 

•  Spectral  Shape  Parameters  (controls  whether  spectrum  is  a  JONSWAP,  Sin¬ 
gle  Peaked  TMA  or  severed  other  shapes) 

•  Length  of  Time  Series 

•  Seed  number  for  randomization  of  linear  wave  superposition 

The  program  creates  a  spectrum  of  waves  based  upon  linear  superposition  of  128 
components  or  spectral  lines.  Each  component  represents  a  certain  wave  height  for 
a  certain  frequency.  The  phase  relationship  between  different  linear  components 
is  random.  The  theory  relating  the  wave  mechanics  to  the  board  motion  is  based 
on  linear  wave  theory  (Dean  and  Dalrymple  1984). 

The  TMA  spectral  shape  (Bouws,  et  al.  1985,  Hughes,  1984)  is  a  depth  depen¬ 
dent  form  of  the  JONSWAP  spectral  shape.  The  JONSWAP  spectrum  is  based  on 
empirical  studies  of  wind  wave  spectra  in  the  North  Sea.  The  TMA  spectrum  uses 
over  2800  wind  wave  spectra  at  various  depths  and  locations  to  empirically  define 
the  JONSWAP  parameters  for  wind  seas  that  have  reached  growth  equilibrium  in 
shallow  water.  The  TMA  spectrum  is  single  peaked. 
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THE  ROOT-MEAN-SQUARE  VELOCITY 


3.1  Theory 

The  root-mean-square  water  particle  velocity  under  irregular  waves  has  been 
derived  for  horizontal  velocities  by  Hughes  (1992).  The  outline  of  his  derivation  is 
reproduced  below. 


The  sea  surface  elevation,  r](x,t),  of  a  linear,  partially  reflected,  unidirectional 
wave  on  a  flat  bottom  can  be  expressed  as  (Dean  and  Dalrymple  1984) 

t/(x,  t )  =  a  cos (fex  —  at)  +  aKr  cos (fex  +  at  +  9)  (1) 


where 


a  =  amplitude  of  the  incident  wave, 
k  =  wave  number,  (=  2 tt/T,  where  L  is  wave  length) 
x  =  horizontal  position, 

a  =  angular  wave  frequency  (=  2ir/T,  where  T  is  wave  period), 
t  =  time, 

Kr  =  reflection  coefficient  (=  ardent /“reflected)*  and 
9  —  phase  shift  due  to  reflection. 


The  first  term  in  Equation  1  represents  the  incident  wave  ,  and  the  second  term 
the  partially  reflected  wave.  For  convenience  the  initial  phase  angle  of  the  incident 
component  is  taken  as  zero  at  t  =  0  and  x  =  0. 

The  fluid  particle  velocities  under  such  a  water  surface  can  be  easily  derived 
from  velocity  potential  theory  as 

aqk  cosh  k  (h  +  z)  , 

u(x,z,t)  =  - — — r - [Kr  cos  (kx  —  at  +  9)  —  cos  (kx  +  at)]  (2) 

a  cosh  fen 

aqk  sinh  k  (h  +  z) 

w(x,z,t)  =  - — — — — — - [Kr  sin  (kx  —  at  4-  9)  —  sin  (fex  +  at)]  (3) 

a  cosh  fen 

where 

u  =  horizontal  component  of  velocity, 
w  =  vertical  component  of  velocity, 
g  =  gravitational  acceleration, 
h  =  water  depth,  and 


24 


z  =  vertical  coordinate,  positive  upward,  with  z  =  0  at  the  still  water  level, 
and  z  —  — /z  at  the  bottom. 


By  linearly  superimposing  many  components  and  making  use  of  the  autocor¬ 
relation  function,  it  is  possible  to  derive  expressions  for  urma  and  wTTni ,  the  root- 
mean-squared  (RMS)  velocity  (Hughes  1992).  (In  this  thesis,  lower  case  letters  u 
and  w  represent  statistics  related  to  the  entire  time  series.  Upper  case  letters  U 
and  W  represent  statistics  related  to  the  peak  velocities,  for  example  UTma  and 
fUmean  ■  ) 


U 


2 

rma 


W 


2 

rma 


cosh2  kn(h  +  z) 
cosh2  knh 

sinh2  kn  (h  +  z) 
cosh2  knh 


(l-2 KTn  cos(2knx  +  en)  +  K2Tn) 
(l  +  2 Krn  cos  (2 knx  +  6n)  +  K2n) 


(4) 

(5) 


The  rms  velocity  is  a  useful  and  important  parameter  in  studies  of  turbulence 
and  sediment  transport.  The  parameters  uTm,  and  wrm,  are  simply  the  statistical 
variances  of  a  velocity  time  series. 

Equations  1,2  and  3  are  given  in  slightly  different  forms  in  the  paper  by  Hughes 
(1992).  The  differences  are  a  result  of  the  choice  of  coordinate  orientation,  that  is, 
whether  the  x-axis  points  to  the  right  or  the  left.  This  difference  drops  out  in  the 
derivation  of  Equations  4  and  5.  Thus  Equations  4  and  5  are  applicable  regardless 
of- the  choice  of  coordinate  axis  orientation.  The  difference  may  become  important 
in  studies  of  phase  angle.  However,  in  this  thesis  it  will  be  shown  that  a  random 
phase  angle  works  as  well  as  a  measured  phase  angle  in  applying  Equations  4  and 
5,  thereby  making  the  discrepency  irrelevant. 

Equations  4  and  5  were  applied  in  this  thesis  by  specifying  the  discrete  wavenum¬ 
bers  (fcn),  discrete  angular  wave  frequencies  (<rn),  and  discrete  amplitudes  (an)  of 
the  measured  incident  wave  spectrum,  along  with  water  depth  ( h )  and  spatial 
position  (X,Z).  The  wave  number  for  a  given  water  depth  and  wave  period  can 
be  calculated  exactly  from  the  linear  wave  dispersion  relation.  However,  Nielsen’s 
(1982)  approximate  equation  solves  for  kn  explicitly  and  is  accurate  to  within  one 
percent  in  the  range  0  <  h/L0  <  0.3  (L0  is  the  deepwater  wavelength  defined  as 
L0  =  gT2 /2n).  Nielsen’s  equation  is  applied  in  the  computer  program  that  solves 
Equations  4  and  5. 
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3.2  Experimental  Results 


Figures  7  to  10  compare  the  velocity  statistics  uTm,  and  wrm,  measured  by  the 
LDV  with  urm,  and  wrm,  estimated  using  the  measured  incident  wave  spectrum  in 
Equations  4  and  5.  Plots  were  made  of  the  ratio  of  measured/predicted  uTTn,  and 
wrmt  vs.  the  following  parameters  of  the  measured  incident  wave  spectrum 

•  Tp , 

•  Hmoi 

•  h,  water  depth, 

•  z/h, 

•  Hmo  /  h, 

•  h/gTp  ,  relative  depth, 

•  Hmo/gTp, 

•  Ursell  Number  =  L20Hmo/h3 , 

•  zHmo/h2. 


where  L0  is  the  deepwater  wavelength  associated  with  the  peak  spectral  period, 
T 

■ip- 


Of  particular  interest  is  the  influence  of  nonlinear  wave  conditions.  Nonlinear- 
ities  may  affect  the  ability  of  Equations  4  and  5  to  predict  uTm,  and  wrms-  A  goal 
of  this  thesis  was  to  develop  empirical  corrections  to  Equations  4  and  5  to  account 
for  any  nonlinearities.  However,  no  clear  trends  are  evident  in  the  plots  shown 
in  Figures  7  to  10  (additional  plots  are  in  Appendix  I,  Figures  37  to  50).  Two 
dimensionless  parameters  commonly  correlated  with  nonlinearities  are  the  Ursell 
number  and  relative  depth.  Figure  7  graphs  measured/predicted  urm»  vs.  relative 
depth.  Hughes  (1992)  produces  the  same  plot  for  the  case  of  perfect  reflection 
from  a  vertical  wall  (rather  than  variable  reflection  from  a  spending  beach  in  this 
study).  Hughes  discovered  a  trend  whereby  measured/predicted  uTrn,  is  under- 
predicted  at  low  values  of  relative  depth  and  is  overpredicted  at  higher  values. 
No  such  trend  is  evident  in  Figure  7.  If  anything,  there  may  be  a  slight  decrease 
for  higher  values  of  relative  depth.  The  author  has  no  physical  explanation  as  to 
why  there  would  be  a  trend  in  the  case  of  a  perfectly  reflecting  vertical  wall  and 
no  trend  for  small  reflection  from  a  spending  beach.  Future  investigations  may 
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urms(measured)/urms(predicted) 
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Figure  8:  urm>  Measured/Predicted  vs.  Ursell  No.  ( HmoL\jh 3) 


wrms(measured)/wrms(predicted) 


Figure  9:  wrmt  Measured/Predicted  vs.  Relative  Depth 


want  to  look  at  intermediate  cases  of  reflection  coefficients  between  that  of  perfect 
reflection  ( Kr  =  1.0)  and  small  reflection  ( Kr  <  20%)  in  order  to  quantify  any 
transition  between  the  trends  discovered  by  Hughes  (1992)  and  the  results  in  this 
thesis. 

The  only  clear  observation  is  that  Equations  4  and  5  slightly  underpredict  the 
measured  urm,  and  wrm,  (except  for  wrmi  near  bottom).  The  reason  for  predicted 
urm,  and  wTTns  being  lower  may  be  due  to  the  effect  of  reflected  waves  “cancelling 
out”  through  linear  superposition  the  incident  wave  velocities,  thereby  lowering 
the  actual  u,  w  and  by  extension  urmt  and  wTmB. 

Figures  7  to  10  show  a  slight  difference  between  points  near  bottom  and  those 
at  mid-depth.  Equations  4  and  5  underpredict  the  magnitude  of  urm ,  and  wrrn, 
more  at  mid-depth  than  at  the  bottom.  This  effect  agrees  with  the  well  known 
observation  that  linear  wave  theory  becomes  an  increasingly  poor  predictor  of 
wave  properties  as  the  wave  crest  is  approached.  The  author  concludes  that  this 
nonlinear  effect  exists  throughout  the  water  column.  However,  the  magnitude 
of  the  difference  becomes  much  less  as  the  bottom  is  approached.  One  would 
intuitively  expect  some  effect  due  to  nonlinearities.  The  differences  for  water 
particle  velocities  can  be  expected  to  parallel  the  differences  for  water  surface 
profile.  However,  the  shape  of  Cnoidal  and  Stokes  profiles  is  such  that  the  variance 
of  the  surface  elevations  rjrm,  is  not  greatly  different  from  linear  theory,  even  though 
individual  water  elevation  points  differ.  One  would  expect  the  same  thing  for 
velocity  measurements  and  this  is  confirmed  by  this  study. 

Plots  which  apply  the  Ursell  number  (Figures  8  and  10)  show  similar  scatter 
compared  to  plots  with  relative  depth  as  the  dependent  variable.  The  lack  of 
any  strong  trends  and  the  inconsistencies  described  in  the  abo/e  paragraph  leads 
the  author  to  conclude  that  development  of  an  empirical  correction  to  Equations 
4  and  5  based  on  the  existing  data  would  not  be  justified.  It  is  concluded  that 
Equations  4  and  5  provide  estimates  correct  to  within  10%  based  on  examination 
of  Figures  7  to  10,  except  for  near  bottom  vertical  velocities.  For  near  bottom 
vertical  velocities,  ( z/h  >  0.85),  the  author  recommends  multiplying  Equation  5 
by  1.5  to  estimate  wrm,  to  within  75%. 

An  important  aspect  making  it  difficult  to  quantify  errors  due  to  nonlinearities 
is  that  the  incident  wave  spectrum  is  being  calculated  by  the  Goda  method,  which 
is  itself  based  on  linear  wave  theory.  Goda  and  Suzuki  (1976)  do  not  rigorously 
quantify  the  errors  in  their  method.  Repeatability  types  of  tests  in  Goda  and 
Suzuki  s  experiments  showed  “dispersion”  of  estimated  incident  amplitudes  about 
the  mean  on  the  order  of  ±5%.  Plots  of  these  tests  demonstrated  random  errors 
with  an  approximately  Gaussian  distribution. 
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Also  of  interest  is  that  Equations  4  and  5  require  specifying  horizontal  position, 
x,  and  phase  angle,  0.  Without  a  reflecting  structure,  KTn  =  0  and  x  and  6n  drop 
out.  For  a  perfectly  reflecting  wall  (and  therefore  presumably  vertical),  6n  drops 
out  and  x  is  the  horizontal  distance  from  the  wall.  However,  for  a  sloping  spending 
beach  with  variable  KTn,  specifying  x  and  0n  is  necessary.  For  such  a  case,  0n  is 
not  always  readily  available  and  x  is  not  precisely  defined.  Figures  11  and  12  plot 
predicted  values  of  uTm,  from  equation  4  vs.  distance  from  toe  of  slope  (values 
of  x).  Note  that  uTm,  does  not  equal  zero  at  x  =  0,  although  a  cursory  look 
at  the  graphs  may  give  that  impression.  Two  lines  are  plotted,  one  calculated 
using  the  measured  phase  angle  and  one  with  a  random  phase  angle.  The  plots 
show  some  interesting  features  such  as  an  envelope  effect  for  the  amplitudes  of 
the  varying  urmt.  Also  evident  is  a  wavelength  on  the  order  of  5.0  m  combined 
with  shorter  wavelengths  on  the  order  of  0.2  m.  The  greatest  variability  in  uTTna  is 
within  one  meter  (one  meter  in  the  wave  tank)  of  the  toe  of  slope.  However,  the 
range  of  variability  does  not  seem  to  decrease  past  one  meter.  These  features  are 
interesting  but  also  peripheral  to  this  thesis.  The  relevant  points  are  as  follows: 


•  The  range  of  variability  of  urmt  and  wTm,  with  varying  x  is  less  than  0.05%, 
for  horizontal  locations  more  than  1  meter  from  reflecting  structure  (in  wave 
tank). 

•  There  is  only  a  minor  difference  between  using  a  random  phase  angle,  (0n), 
and  a  measured  phase  angle  in  applying  Equations  4  and  5.  A  measured 
phase  angle  results  from  the  Fourier  analysis  of  the  wave  signal.  A  random 
phase  angle  is  a  random  value  between  0°  and  90°  applied  in  Equation  4  and 

5. 

•  The  author  recommends  simply  using  the  distance  from  the  toe  of  slope 
for  the  horizontal  position,  x.  For  a  structure  with  a  reflection  coefficient 
above  20%,  this  may  not  be  appropriate.  For  this  case,  variability  of  uTTn, 
with  x  may  become  important,  as  discussed  by  Hughes  (1992).  The  reason 
for  such  variability  can  be  intuitively  understood  by  considering  a  regular 
wave  reflecting  from  a  vertical  wall.  The  clapotis  that  is  established  has 
a  region  half  a  wavelength  from  the  wall  where  the  horizontal  velocity  is 
theoretically  always  zero.  One  quarter  of  a  wavelength  from  the  wall  the 
horizontal  velocity  is  a  maximum.  Thus  a  narrow  banded  spectrum  will 
approach  the  behaviour  of  such  a  regular  wave  resulting  in  a  non-uniform 
distribution  of  uTms  with  x.  A  non-uniform  distribution  of  wave  height  will 
also  occur. 

•  The  greatest  variability  in  predicted  urmt  is  adjacent  to  the  reflecting  struc¬ 
ture. 
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urms,  cm/sec 
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Figure  11:  Predicted  urm.  Velocity  vs.  X-coordinate  (Time  Series  No.  9). 
h  =15  cm,  Z  =25  cm,  Hmo  =4.2  cm,  Tp  =1.0  sec 


urms,  cm/sec 


34 


Figure  12:  Predicted  uTm ,  Velocity  vs.  X-coordinate  (Time  Series  No.  69). 
h  =  15  cm,  Z  =7.5  cm,  Hmo  =3.2  cm,  Tv  =1.0  sec 


Figure  13  plots  the  %  range  of  variation  in  predicted  urmt  (as  it  varies  with  x) 
vs.  relative  depth  for  every  time  series.  Additional  plots  are  in  Appendix  I,  Figures 
51  to  53.  No  trends  are  observed  with  respect  to  the  relative  depth  (h/gT£).  It 
is  evident  that  the  range  of  variability  of  urma  and  xuTrnt  with  respect  to  x  is  very 
small  for  the  case  of  Kr  <  20%. 


3.3  Laboratory  Effects 

Figures  14  to  17  plot  the  reflection  coefficient,  Kr,^k-  various  parameters. 
Figure  14  presents  the  result  of  regular  waves  that  wereMmn  before  the  irregu¬ 
lar  wave  time  series.  These  18  regular  wave  runs  were  made  in  order  to  verify 
some  of  the  computer  programs  and  to  optimize  the  design  of  the  wave  absorb¬ 
ing  (spending)  beach.  The  final  configuration  of  the  spending-beach  used  three 
layers  of  “horsehair”  wave  absorbing  mats  (approximately  2.5  cm  per  mat)  on 
top  of  the  T.30  concrete  slope.  Low  Kr  values  were  desired  for  the  experiments. 
The  relatively  high  Kr  values  reported  here  are  in  some  ways  disappointing.  It 
is  interesting  that  increased  layers  of  horsehair  mats  slightly  increased  Kr  (Figure 
14).  The  author  suspects  that  vertical  wave  screens  and/or  a  gravel  surface  on  the 
sloping  beach  would  have  further  reduced  reflection. 

Figure  14  depicts  a  trend  of  increasing  reflection  for  increasing  ratio  of  wave 
height  to  water  depth.  This  trend  is  even  more  pronounced  for  the  case  of  ir¬ 
regular  waves  (Figure  15).  This  trend  may  be  partly  a  laboratory  effect  caused 
by  the  method  of  measuring  reflection.  Wave  height/water  depth  is  also  directly 
correlated  with  incidents  of  wave  breaking  in  the  wave  tank.  The  random  su¬ 
perposition  of  wave  components  results  in  instances  where  individual  waves  are 
excessively  high  or  steep  resulting  in  wave  breaking.  This  type  of  wave  breaking 
can  occur  anywhere  in  the  wave  tank.  The  breaking  of  the  wave  results  in  sec¬ 
ondary  waves  propagating  in  both  directions  in  the  wave  channel.  Thus,  wave 
energy  that  was  travelling  solely  in  a  direction  incident  to  the  beach  is,  through 
the  breaking  process,  sent  in  both  the  incident  and  reflected  directions  thereby 
increasing  the  measured  value  of  Kr.  Higher  wave  heights  will  result  in  a  greater 
number  of  instances  of  breaking  waves.  Since  most  incident  waves  are  higher  than 
reflected  waves,  the  breaking  process  should  result  in  a  net  shift  of  energy  from 
incident  waves  to  reflected  waves. 

Such  a  laboratory  effect  could  possibly  be  minimized  by  locating  the  Goda 
array  (water  elevation  gauges)  closer  to  the  structure  (spending-beach).  Note, 
however,  that  Goda  and  Suzuki  (1976)  advise  locating  the  wave  gauges  used  for 
wave  resolution  at  a  distance  of  more  than  one  wavelength  from  both  the  test 
structure  and  wave  board.  For  this  reason  the  wave  gauges  used  were  located 
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Figure  13:  %-Range  of  urmi  Predicted  as  a  Function  of  X,  (Random  Phase  Angle) 
vs.  Relative  Depth 
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Figure  14.  Reflection  Coefficient  vs.  Hjh ,  Monochromatic  Waves 
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Figure  15:  Reflection  Coefficient  vs.  Hmo/h ,  Irregular  Waves 
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Figure  16:  Reflection  Coefficient  vs.  Relative  Depth,  Irregular  Waves 
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Figure  17:  Reflection  Coefficient  vs.  Peak  Spectral  Period,  Irregular  Waves 


midway  between  the  wave  board  and  toe  of  slope. 


Working  against  this  hypothesis  of  distortion  of  Kr  by  wave  breaking  is  that 
the  regular  wave  runs  did  not  have  any  observed  wave  breaking  between  the  Goda 
array  and  the  spending  beach  for  the  case  of  H/h  <  0.4 .  Despite  this,  a  trend  is 
still  apparent  (Figure  14). 

It  is  fortunate  that  h/gTp  is  not  closely  correlated  with  Kr  (Figure  16)  since 
this  dimensionless  parameter  is  applied  in  empirical  formulae  later  in  this  thesis. 
Also  of  interest  is  that  Tp  is  not  correlated  with  KT  (Figure  17).  This  is  unusual, 
one  would  expect  a  strong  positive  correlation  between  Tp  and  KT  since  higher 
frequency  waves  are  more  easily  dissipated  by  the  horsehair  wave  absorbing  mats. 

The  plots  show  experiment  repeatability  by  differentiating  between  mid-depth 
and  near  bottom  LDV  measuring  points.  (Indicated  by  labeling  the  data  points 
either  Run  No.  1  or  Run  No.  2.)  Although  this  distinction  is  not  relevant  to  the 
measurement  of  Kr ,  values  of  Kr  are  nearly  identical  for  the  theoretically  identical 
wave  conditions  thereby  showing  good  experiment  repeatability. 

Experiment  repeatability  is  also  indicated  by  time  series  Numbers  340,  581, 
712,  and  880.  These  four  runs  were  the  same  as  previous  runs.  The  %  differences 
for  certain  measured  variables  for  these  four  repeated  runs  are  listed  in  Table  3. 
(U  and  W  refer  to  peak  velocities  defined  by  a  down-crossing  analysis,  rms  and 
rmq  refer  to  root-mean-square  and  root-mean-quad  velocities,  (  +  )  and  (-)  refer 
to  the  positive  and  negative  directions  respectively.  These  parameters  are  applied 
later  in  this  thesis.) 

Spatial  homogeneity  in  the  X-direction  describes  the  quality  whereby  exper¬ 
imental  results  are  independent  of  horizontal  distance  from  the  spending  beach. 
Time  series  numbers  341,  342,  582,  583,  713,  and  714  are  spatial  homogeneity 
tests.  Wave  reflection  from  a  structure  will  cause  some  spatial  variability  in  mea¬ 
sured  velocities  and  wave  height  statistics  for  the  reasons  explained  in  Section  3.1. 
However,  no  consistent  modulation  of  the  predicted  value  of  uTTnt  is  evident  in 
Figures  11  to  12.  For  low  values  of  Kr,  say  less  than  20%,  it  can  be  concluded  that 
the  measured  statistics  will  be  nearly  uniform  despite  varying  X.  The  statistics  of 
the  spatial  homogeneity  tests  support  this  conclusion.  The  observed  variability  in 
the  measured  statistics  for  varying  horizontal  location  (X)  is  comparable  to  that 
for  the  repeatability  tests  in  Table  3. 

A  comparison  of  the  number  of  zero-crossing  “waves”  counted  separately  from 
the  particle  velocity  time  series  and  the  water  elevation  time  series  is  depicted  in 
Figuie  18.  Additional  plots  are  in  Appendix  I,  Figures  54  and  55.  The  number  of 
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Table  3:  Experiment  Repeatability 


Measured 

%  Difference  -  Run  No. 

Variable 

34  k  340 

58  k  581 

71  k  712 

88  k  880 

T 

xv 

0.5% 

0.5% 

0.5% 

0.5% 

Hmo 

0.2% 

0.2% 

0.5% 

0.2% 

Kr 

0.1% 

0.6% 

0.7% 

0.7% 

^rmi 

0.1% 

0.4% 

0.2% 

0.2% 

Wrms 

0.3% 

1.6% 

1.9% 

0.4% 

Urrn.(  +  ) 

0.2% 

3.0% 

0.2% 

0.7% 

Urm.(~) 

0.3% 

0.3% 

0.2% 

0.2% 

wrm,(+) 

1.5% 

8.0% 

23% 

50% 

1.0% 

10% 

1.6% 

2.9% 

Urmqi  +  ) 

0.3% 

9.5% 

0.1% 

0.6% 

Urmq(-) 

0.4% 

1.4% 

6.4% 

0.2% 

wrmq(+) 

3.0% 

7.5% 

3.1% 

8.3% 

h 

30  cm 

20  cm 

15  cm 

15  cm 

z 

mid-depth 

mid-depth 

mid-depth 

near  bottom 

H mo 

8.25  cm 

7.51  cm 

3.33  cm 

6.15  cm 

T 

xv 

1.97  sec 

2.87  sec 

1.54  sec 

3.01  see 

waves  has  been  prorated  to  account  for  the  slightly  different  record  lengths.  Each 
wave  is  by  two  consecutive  downcrossings  of  the  mean.  The  number  of  velocity 
waves  measured  by  the  LDV  and  the  number  of  water  elevation  waves  measured 
by  the  wave  gauges  should  theoretically  be  the  same.  For  horizontal  velocities,  the 
number  of  waves  counted  was  approximately  20%  less  than  the  number  of  surface 
waves.  For  vertical  velocities,  some  interesting  trends  are  evident  with  respect  to 
increasing  nonlinearities  (discussed  further  below). 

The  velocity  waves  were  counted  after  the  time  series  had  undergone  a  band¬ 
pass  filter  of  0.20-5.00  seconds  (Figure  6).  This  filtering  removed  high  frequency 
oscillations  in  the  velocity  time  series  caused  by  turbulence  and  imperfections  in  the 
electronic  data  acquisition  systems.  Low  frequency  oscillations  were  also  removed. 
The  effect  of  filtering  is  to  decrease  the  number  of  “velocity”  waves  counted.  Note 
that  for  turbulence  studies  it  may  be  desirable  to  retain  the  effects  of  turbulent 
fluctuations. 

Figure  18  provides  a  useful  check  on  the  quality  of  the  measurements  and  the 
effectiveness  of  filtering.  Large  discrepancies  between  wave  counts  for  velocities 
and  water  elevations  is  a  warning  that  something  may  be  out  of  order.  The  author 
concludes  that  all  experiment  measurements  were  high  quality  based  on  these  wave 
counts  as  well  as  a  visual  examination  of  parts  of  every  velocity  time  series  (Figure 
6).  A  discrepancy  of  more  than  50%  could  be  considered  “poor”  quality  ior  this 
type  of  study.  The  author  has  no  analytical  explanation  for  the  trend  observed  in 
Figure  18b  (increasing  number  of  vertical  velocity  waves  with  decreasing  relative 
depth).  Relative  depth  is  commonly  used  to  indicate  nonlinearities.  The  plot  shows 
an  increase  in  the  number  of  “velocity”  waves  with  increasing  nonlinearities.  One 
can  conclude  from  this  that  up/down-crossing  analyses  have  limited  usefulness  in 
studies  requiring  precise  statistics  or  attempting  to  show  subtle  effects. 
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Figure  18:  Number  of  Waves  Counted  for  Water  Elevation  Time  Series  and  Water 
Particle  Velocity  Time  Series:  (a)  Horizontal  Velocities;  (b)  Vertical  Velocities 
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4  VELOCITY  DISTRIBUTION 


Linear  irregular  wave  theory  assumes  that  the  time  series  of  sea  surface  ele¬ 
vations  is  the  result  of  linear  superposition  of  many  sinusoidal  wave  components, 
each  expressed  by  Equation  1.  That  is,  each  component  has  a  unique  frequency, 
amplitude,  wave  length,  phase  angle  and  reflection  coefficient  and  these  compo¬ 
nents  do  not  interact.  Also  assumed  is  that  sea  surface  elevation,  77 (f),  at  a  fixed 
point  in  space  is  a  random  process  with  a  Gaussian  (Normal)  distribution. 

The  time  series  of  fluid  particle  velocity  also  can  be  assumed  to  be  the  result 
of  linear  superposition  of  many  sinusoidal  components.  Each  component  being 
expressed  by  Equation  2  or  3.  Also,  it  can  be  reasonably  assumed  that  u(t) 
and  w(t)  follow  a  Gaussian  distribution  with  zero  mean  and  standard  deviation 
(urm,  and  wTmt)  given  by  Equations  4  and  5  (Hughes  1992).  The  development 
of  a  probability  distribution  for  particle  velocity  will  parallel  that  for  sea  surface 
elevation. 


4.1  Gram-Charlier  Distribution 


Longuet-Higgins  (1963)  developed  a  probability  distribution  function  applica¬ 
ble  to  non-Gaussian,  nonlinear  ocean  waves  through  a  modification  of  the  normal 
distribution  (O.S.  Rice  (1944)  originally  developed  much  of  the  theories  relating  to 
the  mathematical  description  of  wave  phenomenon).  The  equation  (in  standard¬ 
ized  form,  after  dividing  by  the  standard  deviation  of  x)  is  given  by  Ochi  (1982) 
as  follows: 


/(*)  = 


exp(~r) 

\Z2tt 


1  +  ^7*3  (*)  +  ^n,(x)  +  + 


3! 


4! 


72 


5! 


(6) 


where 
x  = 

Hn  = 

^3  = 

A4  -- 

A  5  = 


nondimensional  parameter  given  by  (x  —  x) / a  where  “x”  is  the  variable 
of  interest  (t/,u,  etc.),  and  <7  is  the  standard  deviation  of  the  variable. 
Hermite  polynomial  of  degree  n, 
skewness  (non-dimensional), 
kurtosis  (non-dimensional), 

higher  order  moment  of  the  nondimensional  distribution. 


This  equation  will  be  referred  to  henceforth  as  the  Gram-Charlier  distribution. 
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Different  probability,  density  functions  are  defined  by  including  different  terms 
of  the  series  in  Equation  6  (Ochi  1982).  The  following  four  distributions  were 
studied: 


Gaussian  distribution:  First  term  of  the  series  in  Equation  6 
Nonlinear  (I):  Up  to  the  third  term  of  the  series  in  Equation  6 
Nonlinear  (II):  Up  to  the  fourth  term  of  the  series  in  Equation  6 
Nonlinear  (III):  All  five  terms  of  the  series  in  Equation  6 


Note  that  the  Hermite  polynomials  used  in  Equation  6  are  not  the  ones  typically 
found  in  math  references  (CRC  Standard  Mathematical  Tables  1987).  Longuet- 
Higgins  (1963)  applies  the  following  version 


Hn 

n  n  (n  —  1)  xn~2  n  (n  —  1)  (n  —  2)  (n  —  3)  x"~4 

1!  2  +  2!  22 

,  (7) 

h3 

=  x3  —  3x, 

(8) 

H< 

=  x4  —  6x2  4-  3, 

(9) 

h5 

=  x5  —  10x3  4-  15x, 

(10) 

H6 

=  x6  -  15x4  +  45x2  -  15. 

(11) 

Application  of  the  Hermite  polynomials  given  in  the  CRC  Tables  does  not  work. 

The  results  of  previous  investigations  using  this  distribution  are  discussed  in 
Section  1.1  of  this  thesis.  Ochi  (1982)  outlines  the  derivation  of  Equation  6,  which 
involves  application  of  the  cumulants-generating  function. 

Note  that  Equation  6  is  close  to  a  similar  equation  using  Gram-Charlier  series 
(Cramer  1946)  defined  as 

/(*)  = 


exp 


(¥) 


y/2^ 


1  I  W  f'-r'l  I  -  3)  rr  /_\  ,  (A6  ~  10A3) 

1  4-  —  Hz\x)  + - — - H4(x)  + - — - H  5(1)  4- 


(12) 


Application  of  Equation  6  to  the  measured  velocity  data  results  in  a  worse  fit 
to  the  data  than  application  of  the  plain  Gaussian  distribution.  Figures  19  and 
20  show  typical  results  for  this  study.  Additional  plots  are  in  Appendix  I,  Figures 
73  to  96.  The  theoretical  probability  distribution  underestimates  the  number  of 
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higher  magnitude  velocity  points  and  overestimates  the  number  of  points  near 
zero.  Figure  4  is  a  sample  hardcopy  output  from  the  LDV  microcomputer.  It 
includes  the  parameters  used  to  calculate  the  theoretical  distributions  as  well  as 
tabulations  of  the  measured  and  theoretical  distributions.  The  tabulated  data  is 
plotted  in  graphs  like  Figures  19  and  20.  Under  no  conditions  does  the  Gram- 
Charlier  distribution  provide  a  better  fit  than  the  Gaussian  distribution. 

Figure  21  plots  %-RMS  vs.  relative  depth  for  each  time  series.  The  value  of  %- 
RMS  error  is  calculated  for  each  time  series  as  follows;  The  difference  between  the 
theoretical  and  measured  probability  density  for  every  non-dimensional  velocity 
point  is  calculated.  The  differences  are  squared  and  then  the  sum  of  these  squares 
is  divided  by  the  number  of  points.  The  square  root  of  this  number  is  the  RMS 
error.  Multiplying  the  RMS  error  by  the  number  of  points  gives  the  %-RMS  error. 
Using  the  %-RMS  error  instead  of  the  plain  RMS  error  is  not  really  necessary 
for  comparing  the  Gram-Charlier  distribution  applied  in  this  thesis.  However,  the 
%-RMS  error  needs  to  be  used  for  the  case  of  an  inconsistent  number  of  points, 
such  as  for  the  Beta-Rayleigh  distribution  later  in  this  thesis. 

The  velocity  time  series  applied  in  this  thesis  were  unfiltered  (raw  data).  Fig¬ 
ures  22  and  23  compare  theoretical  and  measured  probability  distributions  for 
filtered  and  unfiltered  data.  The  differences  are  not  significant,  especially  when 
comparing  the  relative  differences  between  the  Gaussian  (linear)  and  nonlinear 
distributions.  Additional  plots  are  contained  in  Appendix  II. 

One  question  is  whether  the  author  has  implemented  the  Gram-Charlier  distri¬ 
bution  correctly.  Unfortunately,  the  expression  does  not  transform  itself  into  the 
plain  Gaussian  distribution  for  certain  theoretical  values  of  the  parameters  A3,A4, 
and  A5.  This  makes  verification  difficult.  The  fact  that  the  probability  distribution 
sums  to  1.0  (Figure  4)  is  one  check  on  the  implementation  although  not  foolproof. 
Another  check  is  the  results  of  previous  investigations  (Section  1.1).  Ochi  (1982) 
reproduces  the  results  of  one  time  series  made  by  Honda  and  Mitsuyasu  (1976) 
for  wave  accelerations.  The  relative  shape  of  the  probability  distribution  curves  is 
very  similar  to  the  shapes  in  this  thesis. 

Since  the  Gram-Charlier  distribution  fits  poorly  the  measured  velocity  distri¬ 
bution,  no  attempt  has  been  made  to  relate  the  parameters  of  the  distribution  to 
parameters  of  the  wave  field. 
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Figure  19:  Horizontal  Velocity  Probability  Distribution 
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Figure  20:  Vertical  Velocity  Probability  Distribution 
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Figure  21:  %-Root-mean-square  Error,  Gram-Charlier  Distribution  vs.  Hmo/h.  (a) 
Horizontal  Direction;  (b)  Vertical  Direction 
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Figure  23:  The  Effect  of  Varying  Filter  Bandwidths  on  Fitting  the  Gram-Charlier 
Distribution  to  Vertical  Velocities,  Time  Series  No.  37.  (a)  Raw  Data;  (b)  Filtered 
Data  [0.20-5.00]  s;  (c)  Filtered  Data  [0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 


5  PEAK  VELOCITY  DISTRIBUTION 


Predicting  the  magnitude  of  wave  height  as  a  function  of  statistical  parameters 
has  been  a  fundamental  goal  of  much  research  involving  random  seas  (wave  height 
is  defined  as  the  vertical  distance  from  a  wave  crest  to  a  successive  trough).  It  can 
be  derived  theoretically  that  the  statistical  distribution  of  wave  heights  (surface 
waves)  follows  the  Rayleigh  probability  law  if  one  assumes  the  following  (Longuet- 
Higgins  1952): 


•  The  time  series  of  water  surface  elevations  is  the  result  of  the  linear  super¬ 
position  of  many  sinusoidal  components. 

•  At  a  fixed  location,  the  water  surface  elevation  follows  a  Gaussian  distribu¬ 
tion. 

•  The  spectrum  of  water  elevations  is  narrow  banded.  That  is,  the  frequencies 
of  the  sinusoidal  components  are  close  together. 

•  The  wave  peaks  are  statistically  independent. 


The  Rayleigh  distribution  has  a  good  history  of  predicting  wave  heights  even  for 
sea  conditions  well  outside  the  theoretical  limits  listed  above.  This  one-parameter 
distribution  is  defined  as 


P(H)  = 


(13) 


where 

P(H)  =  Probability  density  function  for  wave  heights, 
H  —  Wave  height, 

HTTn,  =  Root-mean-squared  wave  height. 


5.1  The  Beta  Rayleigh  Distribution 


Hughes  and  Borgman  (1987)  applied  the  “Beta- Rayleigh”  distribution  that 
modifies  the  above  equation  as  follows 

0-1 


2a  —  1 


Pbr  (H )  = 


2T  (a  +  0)  H 
T(a)T(0)  Hf 


1  - 


fp  v 

m) 


0  <  H  <  Hh 


(14) 
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where 


Hb  ~  Breaking  wave  height  (the  maximum  wave  height  possible), 
T  (  •  )  =  The  Gamma  function. 


a  = 

JG  (K2  -  Ki) 

K\  -  K2 

(15) 

0  = 

i 

i 

r—H 

(16) 

K\  -  K2 

K\  = 

HL. 

«S 

(17) 

k2  = 

Ht 

(18) 

Hrmq  = 

Root-mean-quad  wave  height 

This  formulation  is  a  three  parameter  distribution  (Hb,  HrmB,  and  Hrmq)  and 
so  should  be  an  improvement  over  the  one  parameter  Rayleigh  distribution.  An¬ 
other  improvement  is  that  the  distribution  is  bounded  by  a  maximum  wave  height, 
whereas  the  Rayleigh  distribution  “tail”  trails  off  into  infinity.  For  deepwater 
waves,  it  is  the  tail  region  that  shows  the  greatest  deviation  between  theory  and 
measured  wave  heights.  Finally,  the  Beta-Rayleigh  distribution  transforms  itself 
into  the  Rayleigh  distribution  when  Hb  approaches  infinity  and  HTTnq  =  \Z2Hlms. 

Comparison  of  measured  wave  heights  with  those  predicted  by  the  Beta-Rayleigh 
and  Rayleigh  distributions  reveals  that  the  new  formulation  is  a  significant  im¬ 
provement  over  the  plain  Rayleigh  distribution.  Also,  empirical  equations  have 
been  developed  relating  the  three  parameters  of  the  distribution  to  the  following 
common  variables  (Hughes  and  Borgman  1987): 

Tp  —  peak  spectral  period, 

h  =  water  depth, 

Hmo  —  zeroth-moment  wave  height,  equal  to  four  times  the  standard  deviation 
of  sea  surface  elevations. 

This  thesis  applies  the  Beta-Rayleigh  distribution  to  particle  velocities  under 
irregular  shallow  water  waves  instead  of  wave  heights.  For  the  time  series  of 
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horizontal  velocities,  the  following  one  for  one  substitutions  are  made  resulting  in 
Equation  23. 


u 

for 

V, 

^ rma 

for 

Vrm, » 

£7/2 

for 

H, 

UTm,/ 2 

for 

HTma  i 

UTmq/2 

for 

II 

**Tmq  i 

Umax/  2 

for 

II 

max  • 

Similar  substitutions  are  made  for  vertical  velocities,  ( w  and  W).  The  peak 
velocity  is  divided  by  2  to  account  for  the  fact  that  peak  particle  velocities  are 
measured  from  the  mean  velocity  level  to  the  maximum  positive  or  negative  ve¬ 
locity  (crest  or  trough  to  mean  level).  Wave  heights,  H,  are  measured  from  the 
crest  to  the  trough. 

For  the  Rayleigh  distribution,  the  following  theoretical  relationships  hold 
(Htghes  1992,  Hughes  and  Borgman  1987) 


Urma 

=  V2.Ur  m. 

(19) 

Urmq 

=  2\/r2u2  . 

v  rm# 

(20) 

The  following  parameters  can  be  determined  from  measurements  using  the 
relationships 

Urm 5  =  (n£U?)  (21) 

( i  n  \  ^ 

v™,  =  (22) 

U  =  peak  velocity  (either  positive  or  negative  direction) 

Again,  similar  equations  apply  for  vertical  velocities. 

Note  that  in  studying  peak  velocity,  one  usually  distinguishes  between  velocity 
in  the  positive  and  negative  direction,  unlike  wave  hfght  that  is  defined  by  the 
total  distance  between  peak  and  trough. 


55 


The  Beta-Rayleigh  probability  distribution  function  applicable  to  peak  particle 
velocities  is  therefore 


Pbr(U)  = 


2T  {a.  +  fi)  U2a 

T(a)rw)  u* 


2  \ 


u 2 

max 


where 


KX(K2-KX) 
a  Kl  -  K2 

(1  -KiHKt-Ki) 

P  ~  K\-  K2 


K  UL. 

Kl  =  uiZ 


K  ULq 
K*  ~  uiZ 


with  similar  equations  applying  for  vertical  velocities. 


5.2  Experimental  Results 


Practical  application  of  the  Beta-Rayleigh  distribution  can  be  done  only  if 
the  three  parameters  of  the  distribution  (Urpax,  Urm,,  and  Urmq)  are  expressed  in 
terms  of  common  parameters  of  the  wave  field.  For  example,  Hughes  and  Borgman 


(1987)  developed  empirical  expressions  relating  Hb,  Hr 
and  h  for  the  case  of  shallow  water  irregular  waves. 


,,  and  Hr mq  with  H„ 


This  investigation  uses  laboratory  measurements  to  develop  empirical  relation¬ 
ships  for  shallow  water  irregular  wave  peak  particle  velocities.  It  will  be  shown 
that  the  empirical  'fit”  is  very  good.  Unfortunately,  the  Beta-Rayleigh  distribu¬ 
tion  provides  only  marginal  improvement  over  the  Rayleigh  distribution. 


Figure  24  plots  measured  and  theoretical  distributions  for  one  of  the  time 
series  (additional  time  series  plots  are  in  Appendix  I,  Figures  73  to  96).  Figure  24 
clearly  shows  the  Beta-Rayleigh  distribution  being  “pulled”  in  ti.?  direction  of  the 
peak  of  the  measured  velocity  distribution.  A  hardcopy  output  from  the  computer 
program  that  calculated  the  statistics  is  in  Figure  4.  The  measured  and  theoretical 
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ure  24:  Velocity  Peaks  Probability  Distribution,  (a)  Horizontal;  (b)  Vertical 


Table  4:  Beta-Rayleigh  Distribution,  Mean  %-RMS  Error 


z=near  bottom 

z=mid-depth 

Direction 

Rayleigh 

Beta-R 

Diff. 

Rayleigh 

Beta-R 

Diff. 

Positive  U 

40.5% 

39.3% 

1.2% 

43.4% 

42.0% 

1.4% 

Negative  U 

40.8% 

36.8% 

2.0% 

44.0% 

38.2% 

5.8% 

Positive  W 

56.1% 

44.4% 

11.7% 

42.8% 

37.3% 

5.5% 

Negative  W 

64.9% 

44.8% 

20.1% 

46.8% 

39.7% 

7.1% 

distributions  tabulated  in  Figure  4  are  an  example  of  what  was  plotted  for  every 
time  series. 

Figures  25  and  26  plot  the  %-RMS  error  for  every  time  series  as  a  function 
of  h/gTp.  Additional  plots  are  contained  in  Appendix  I,  Figures  56  to  58  (The 
use  of  %-RMS  error  as  a  measure  of  “goodness-of-fit”  is  discussed  in  section  3.1). 
No  significant  trends  are  evident  in  the  plots  that  would  define  when  the  Beta- 
Rayleigh  distribution  is  an  improvement  over  the  Rayleigh  distribution.  However, 
on  average  the  Beta-Rayleigh  distribution  is  better  than  the  Rayleigh  distribution 
(Table  4).  Table  4  lists  the  %-RMS  error  for  both  distributions  and  for  different 
velocity  directions.  The  relative  improvement  is  greater  for  near  bottom  vertical 
velocity  than  for  near  bottom  horizontal  velocity. 

Filtering  of  the  velocity  time  series  was  done  before  calculation  of  the  theoreti¬ 
cal  and  measured  distributions.  A  band-pass  filter  in  the  range  [0.20-5.00]  seconds 
was  selected  based  on  comparison  with  the  results  of  other  filter  widths  (Figures 
27,  28  and  in  Appendix  I,  Figures  59  to  64).  Figures  27  and  28  show  that  this 
filter  range  achieves  the  best  agreement  between  the  theoretical  and  measured 
distributions.  However,  a  filter  range  of  [0.20-5.00]  s  was  not  always  optimal.  Fil¬ 
tering  can  make  a  dramatic  difference  in  data  with  high  frequency  turbulence  and 
micro-eddies.  Filtering  also  can  improve  a  “noisy”  signal  caused  by  a  poorly  tuned 
LDV  apparatus  or  other  data  acquisition  imperfections.  Fortunately,  about  half 
the  time  senes  looked  almost  as  good  as  Figure  6b  prior  to  any  filtering.  Applying 
the  filter  to  that  relatively  good  signal  will  not  change  the  down-crossing  analysis 
results  significantly.  Therefore,  in  order  to  be  consistent,  a  filter  of  [0.20-5.00]  s 
was  applied  to  all  time  ..eries  regardless  of  whether  it  was  needed.  Note  that  for 
turbulence  studies  it  may  be  desirable  to  retain  the  effects  of  turbulent  flow  since 
turbulence  is  a  real  phenomena.  However,  this  thesis  assumes  irrotational,  inviscid 
flow  for  which  turbulence  does  not  theoretically  exist. 

Selection  of  bar  width  in  plotting  the  measured  distributions  was  a  challenge. 
Ideally,  one  desires  a  constant  bar  width  for  ease  of  comparison  between  time  series. 
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%-rms  error 


Figure  25:  %-Root-mean-square  Error,  Beta- Rayleigh  Distribution  and  Rayleigh 
Distribution,  Negative  Horizontal  Direction,  Mid-depth 


Figure  26:  %-Root-mean-square  Error,  Beta- Rayleigh  Distribution  and  Rayleigh 
Distribution,  Negative  Horizontal  Direction,  Near  Bottom 
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Figure  27:  The  Effect  of  Varying  Filter  Bandwidths  on  the  Probability  Distribution 
of  Horizontal  Velocity  Peaks,  (a)  Raw  Data;  (b)  Filtered  Data  [0.2-5.00]  s;  (c) 
Filtered  Data  [0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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Figure  28:  The  Effect  of  Varying  Filter  Bandwidths  on  the  Probability  Distribution 
of  Vertical  Velocity  Peaks:  (a)  Raw  Data;  (b)  Filtered  Data  [0.20-5.00]  s;  (c) 
Filtered  Data  [0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 


The  plots  of  the  velocity  distribution  for  the  entire  time  series  have  constant  bar- 
width  (Figures  19  and  20)  because  the  large  number  of  data  points  collected  (8192 
per  time  series)  provide  “stable”  statistics  that  permit  considerable  flexibility  in 
selecting  bar  width.  However,  the  number  of  velocity  peak  data  points  measured 
per  time  series  was  100  to  300.  This  relatively  low  number  of  points  is  adequate  for 
calculating  stable  statistics  but  barely  adequate  for  graphing  the  measured  data. 
Hence,  a  significant  portion  of  the  %-RMS  error  calculated  for  every  time  series 
is  attributed  to  these  errors  in  the  “discretization”  of  the  theoretically  smooth 
probability  distribution. 


5.3  Parameterization 


An  interesting  way  of  comparing  the  Bet  a- Rayleigh  to  the  Rayleigh  distribu¬ 
tion  is  found  by  plotting  Urms/uTm»  and  UrTnq/ulma  vs.  various  non-dimensional 
parameters  (Figures  29  and  30).  These  theoretical  ratios,  assuming  the  Rayleigh 
distribution  (and  its  underlying  assumptions)  holds,  are  constants  given  by  Equa¬ 
tions  19  and  20.  For  the  case  of  horizontal  velocity  there  is  a  clear  divergence  from 
these  theoretical  ratios  for  increasing  nonlinearities  (indicated  here  by  increasing 
Ursell  number).  One  can  conclude  that  horizontal  velocity  peaks  are  greater  in  the 
direction  of  wave  propagation  (in  the  positive  direction)  than  in  the  direction  op¬ 
posite  of  wave  travel  (in  the  negative  direction).  No  such  difference  is  observed  for 
the  case  of  vertical  velocity  peaks.  These  observations  agree  with  nonlinear  second 
order  Stokes  theory  that  predicts  higher  horizontal  velocities  in  the  direction  of 
wave  travel  (which  is  under  the  crest)  than  in  the  opposite  direction  (Dean  and 
Dalrymple  1984).  The  results  do  not  agree  with  previous  investigations  (Section 
1.1).  Note  that  linear  wave  theory  predicts  a  symmetrical  distribution  of  velocity. 
Quantification  of  this  divergence  may  be  worthwhile  but  is  not  directly  relevant 
to  this  thesis. 

Comparison  of  the  Beta- Rayleigh  and  Rayleigh  distributions  is  also  done  by 
plotting  Beta-Rayleigh  distributions  with  the  parameters  (7rmg  and  UTTnq  calculated 
from  Equations  19  and  20  (Figure  31  and  in  Appendix  I,  Figure  65).  The  effect  is 
to  create  a  two  parameter  distribution  based  on  urmt  and  Umax  (wTrnil  and  Wmax  in 
the  vertical  direction).  Clearly,  if  £/max  approaches  infinity  then  this  two  parameter 
distribution  should  transition  into  the  Rayleigh  distribution.  Also,  there  should 
be  no  difference  between  the  positive  direction  and  negative  direction  theoretical 
distributions  since  urmg  and  f7max  apply  to  both  directions.  These  features  are 
evident  in  Figure  31  which  shows  the  Beta-Rayleigh  distribution  nearly  the  same 
as  the  Rayleigh  distribution. 
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Figure  29:  Theoretical  Ratios,  Beta-Rayleigh  Distribution,  Horizontal  Direction, 
(a)  (JrTnt/uTmt ;  (b)  ^  mq/Wrrn, 
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Figure  30:  Theoretical  Ratios,  Beta-Rayleigh  Distribution,  Vertical  Direction,  (a) 

^rmi/^rmn  (b)  Wrmq/WTTnl 
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Figure  31:  Horizontal  Velocity  Peaks  Probability  Distribution,  uTTnt  Only  Param¬ 
eter 


5.3.1  Maximum  Velocity 


The  development  of  an  expression  for  maximum  velocity  is  more  difficult  than 
determination  of  a  breaking  wave  height.  Many  investigations  both  theoretical  and 
empirical  have  been  made  into  the  subject  of  breaking  wave  height.  The  simple 
expression 

Hb  =  0.78A  (28) 

where  h  =  water  depth  is  widely  used  (McGowan  1894). 


The  Shore  Protection  Manual  (1984)  suggests  the  following  expression  for  the 
maximum  horizontal  velocity  under  a  breaking  wave 

Umax  =  Jg(H  +  h)  (29) 


where  (H  +  h)  is  the  distance  measured  front  the  crest  of  the  breaker  to  the  bottom. 


Equation  29  is  based  on  laboratory  measurements  by  Daily  and  Stephan  (1953) 
for  solitary  waves.  The  expression  is  probably  more  applicable  to  wave  action  in 
the  surf  zone  where  a  broken  wave  moves  forward  as  a  bore  of  water.  This  thesis  is 
more  concerned  with  fluid  particle  velocities  under  progressive  waves  where  wave 
breaking  is  governed  more  by  wave  steepness  than  by  limited  depth. 


For  a  shallow  water  regular  wave,  linear  wave  theory  predicts  the  following 
horizontal  and  vertical  velocities  (Shore  Protection  Manual  1984) 

u  =  7T\/?cos^  (30) 

2  V  h 

w  =  ^1  +  ^  sin 9  (31) 

where  shallow  water  is  defined  as 


h  J_ 
L  <  25 


(32) 


6  is  the  phase  of  the  regular  wave  and  z  is  positive  upwards  measun  <1  fr*un  the 
Stillwater  level.  These  expressions  are  derived  from  Equations  2  and  t 


Expressions  for  the  maximum  velocities  can  be  developed  by  a>Mumug  that 
Hmax  —  h  (an  assumption  based  on  Equation  28).  Also  consider  that  the  maximum 
velocity  will  occur  when  cos#  =  1.0.  Finally,  a  reasonable  choice  for  T  in  Equation 
31  is  Tp  Making  these  substitutions  one  arrives  at  the  following  expressions  for 


(33) 


the  maximum  horizontal  and  vertical  velocity 


U„ 


2 


(34) 


A  change  is  made  to  Equation  34  in  applying  it  to  this  investigation.  The 
mean  period  (T)  of  the  velocity  time  series  is  used  instead  of  the  peak  period 
( Tp )  measured  from  the  water  elevation  time  series.  Other  than  this,  Equations 
33  and  34  are  identical  to  those  applied  in  this  investigation.  This  change  was 
made  in  order  to  simplify  the  calculation  procedure.  This  way  all  the  statistics 
of  the  velocity  time  series  could  be  calculated  independent  of  the  water  elevation 
statistics.  The  mean  period  of  velocity  is  generally  less  than  the  peak  period  of 
wave  heights,  at  times  50%  less  than  Tv.  As  will  be  shown  below,  the  effect  on  the 
resulting  theoretical  distributions  is  not  significant. 

Comparison  of  measured  with  theoretical  maximum  velocity  is  made  in  Figure 
32  and  in  Appendix  I,  Figure  66.  A  clear  trend  is  evident  whereby  the  ratio  of 
measured  to  predicted  maximum  velocity  increases  with  Hmo/h.  This  trend  is 
interesting  but  is  not  applied  in  this  thesis.  Figure  33  and  in  Appendix  I,  Figures 
67  and  68  show  that  a  precise  determination  of  Umax  and  Wmax  is  not  necessary 
in  applying  the  Beta-Rayleigh  distribution  to  velocities.  The  difference  in  the 
theoretical  distributions  for  0.5 x  and  2.0 x  the  values  in  Equations  33  and  34  is 
barely  discernible.  A  value  less  than  0.5 x  the  usual  value  of  Lmax  and  Wmax  would 
be  hard  to  justify  in  physical  terms  and  would  be  past  the  limit  of  measured  values 
(Figure  32).  Additional  research  is  probably  worthwhile  on  the  general  subject  of 
“extremal”  velocity  estimates. 

The  assumption  of  “shallow  water”  as  defined  by  Equation  32  deserves  closer 
scrutiny.  Equat' m  33  is  independent  of  z  and  this  is  correct  in  a  physical  sense 
for  shallow  water  linear  waves.  However,  some  components  of  the  wave  spectrum 
probably  do  not  satisfy  Equation  32  and  therefore  some  decrease  in  Umax  with 
increasing  depth  below  water  surface  can  be  expected.  Unfortunately,  the  exper¬ 
imental  program  was  not  intended  to  investigate  fully  the  variation  of  velocities 
with  z  and  therefore  it  is  difficult  to  quantify  the  relationships.  Half  the  measure¬ 
ments  were  taken  at  mid-depth  (z/h  =0.5)  and  half  at  “near  bottom”  (z/h  =0.85 
to  0.94).  As  explained  in  the  above  paragraph,  a  precise  determination  of  Umax  is 
not  important  in  applying  the  Beta-Rayligh  distribution. 

Several  clues  are  available  which  show  that  the  shallow  water  assumption  of 
u  being  independent  of  z  is  reasonable.  First,  Figure  32  plots  the  ratio  of  mea- 
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Umax(measured)/Umax(predicted) 


LI  0.2  0.3  0.4 

Hmo/h 

U,  z  a  mid-depth  '  U,  z  a  near  bottom 


(a)  Horizontal  Direction 


Wmax(measured)/Wmax(predicted) 


0.2  0.3 

Hmo/h 


W,  z  a  mid-depth 


W,  z  a  near  bottom 


b )  Vertical  Direction 


Figure  32:  Maximum  Velocity  Measured  Predicted,  (a)  Horizontal  Direction;  (b) 
Vertical  Direction 
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Figure  33:  Horizontal  Velocity  Peaks  Probability  Distribution  (a)  2x  Usual  Umax 

lb)  0.OX  Usual  Umax 
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sured/predicted  f/max.  No  large  difference  is  evident  between  mid-depth  and  near 
bottom  measurements,  predicted  Umax  is  nearly  the  same  for  both  depths.  One 
would  expect  mid-depth  ratios  to  be  greater  if  u  were  a  function  of  z  (the  same 
wave  conditions  were  run  twice,  one  at  mid-depth  and  one  near  bottom). 

A  second  indication  is  provided  by  Figure  34,  which  plots  the  mean  peak  hor¬ 
izontal  velocity  ( U )  in  the  positive  and  negative  directions  for  every  time  series. 
Again,  one  would  expect  the  mid-depth  horizontal  velocities  to  be  higher  than  the 
near  bottom  velocities  if  the  shallow  water  assumption  is  not  valid.  Here,  the  mid¬ 
depth  U  ts  larger  than  the  near  bottom  U,  but  only  slightly.  A  third  indication  is 
provided  by  Figure  35  (discussed  in  the  next  section).  Again,  there  is  no  practical 
difference  between  near  bottom  and  mid-depth  horizontal  velocity  measurements. 

The  maximum  vertical  velocity,  Wmax ,  predicted  by  Equation  34  is  a  function 
of  z.  However,  the  predicted  values  of  Wmax  are  also  affected  by  the  shallow  water 
approximation  in  linear  theory.  The  effect  is  just  not  as  clear  as  that  for  horizontal 
velocity  where  z  drops  out  completely.  The  discrepancy  due  to  just  the  shallow 
water  approximation  is  difficult  to  quantify  for  w.  This  discussion  of  the  shallow 
water  approximation  is  part  of  a  more  general  concern  regarding  the  range  of 
applicability  of  the  equations  in  this  thesis  (ranges  of  applicability  are  defined  in 
the  next  section). 


5.3.2  Root-Mean-Square  and  Root-Mean-Quad  Peak  Velocity 


Nondimensional  plots  of  the  Root-mean-square  and  Root -mean- quad  peak 
velocities  vs.  relative  depth  are  presented  in  Figures  35  and  36.  Additional  plots 
using  other  nondimensional  parameters  are  in  Appendix  I,  Figures  69  to  72.  Clear 
trends  are  evident  which  lend  themselves  to  empirical  curve  fitting.  Regression 
analysis  yields  the  following  equations 


Urm,Tp 

H  (  +  ) 

“  mo 

=  0.006028  (i)-1  9830  exp 

-0.14391  •  In  (x)2 j 

’  h 

<  0.5 

(35) 

r  t 

c  rms  ± p  ,  v 

IT  '  ' 

**  mo 

=  0.010565  (i)"1  6598  exp 

[-0.09583  -  In (x)21 

i’  h 

<  0.5 

(36) 

r  t 2 

L  rmq  1p 

H 2  ^ 

mo 

=  0.000057  (x)“3  9741  exp 

[-0.29475  -ln(x)2' 

’  h 

<  0.5 

(37) 

IrmqT* 

HI,  ' 

=  0.00  0  2  2  9  (  x )  3  1383  exp 

-0.17244  •  In  (x  )2 

■  ’  h 

<  0.5 

(38) 
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W  T 

Hjjxo 

=  0.178386  (x)-0'4841  exp 

-0.01617  -  In  (x)2' 

,  i=  0.5 

(39) 

Wrm9Tp 

H-mo 

=  0.029270  (x)_0-5147exp 
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h 

(40) 
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h 

x  = 

gTl 

all  valid  only  in  the  range 


0.002  < 


h 

7^ 


<  0.05  . 


The  positive  horizontal  direction  is  defined  as  being  in  the  direction  of  wave 
travel.  The  negative  horizontal  direction  is  opposite  the  direction  of  wave  travel. 
The  positive  vertical  direction  is  upwards. 

Each  equation  is  of  the  form 

y  =  a(x)b  exp  je  •  In  (x)2j  . 

Reasonable  upper  and  lower  limits  to  the  observed  data  are  defined  by  replacing 
"a"  in  Equations  35  through  42  with  the  values  in  Table  5.  Included  in  Table  5 
are  the  correlation  coefficients  for  the  above  equations. 

An  interesting  phenomenon  is  that  for  horizontal  velocities  (Figure  35),  the 
mid-depth  and  near  bottom  measured  data  coincide,  while  the  positive  and  neg¬ 
ative  directions  diverge  for  increasing  nonlinearities.  For  vertical  velocities  (Fig¬ 
ure  36),  the  mid-depth  and  near  bottom  measured  data  are  different  (including 
increased  divergence  with  increasing  nonlinearities)  while  the  positive  and  negative 
directions  coincide.  These  observations  agree  with  linear  wave  theory. 

In  principle,  it  would  be  desirable  to  combine  Equations  39  and  40  as  well 
as  Equations  41  and  42  into  single  equations  that  are  a  function  of  z.  Unfortu¬ 
nately.  the  data  set  only  includes  points  at  mid-depth  and  near  bottom  making 
development  of  such  equations  difficult.  The  author  recommends  applying  linear 
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Table  5:  Root -Mean- Square  and  Root-Mean-Quad  Curve  Fitting 


Equation 

Number 

Correlation 

Coefficient 

a 

upper  limit 

a 

lower  limit 

35 

0.975 

0.007447 

0.004992 

36 

0.987 

0.012916 

0.008651 

37 

0.973 

0.000075 

0.000045 

38 

0.973 

0.000335 

0.000203 

39 

0.942 

0.223130 

0.135335 

40 

0.817 

0.049787 

0.018316 

41 

0.949 

0.082085 

0.040762 

42 

0.832 

0.030197 

0.004087 

interpolation  for  values  of  z/h  between  0.80  and  0.50  (In  principle  the  interpola¬ 
tion  should  be  a  parabolic  function  to  match  linear  wave  theory.  However,  given 
the  scatter  in  the  measured  data,  the  error  in  using  linear  interpolation  should  be 
insignificant.)  For  values  of  z/h  between  0.95  and  0.80,  use  of  Equations  40  and 
42  only  is  recommended.  Application  of  these  equations  above  z/h  =  0.5  or  below 
z/h  =  0.95  is  not  recommended. 


73 


mean  (cm/s) 


CD 

o 

6 


74 


Figure  34:  Mean  Horizontal  Velocity  Peaks 
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Figure  35:  Parameterization  of  Beta-Rayleigh  Distribution,  Horizontal  Direction. 
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Figure  36.  Parameterization  of  Beta- Rayleigh  Distribution,  Vertical  Direction,  (a) 
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6  SUMMARY  AND  CONCLUSIONS 


In  this  study,  irregular  wave  shallow  water  wave  particle  velocities  were  in¬ 
vestigated  with  a  laser  Doppler  velocimeter  (LDV).  Measurements  were  made  in 
a  flat  bottom  2-dimensional  wave  tank  with  a  spending  beach  at  one  end.  Of 
special  interest  is  the  effect  of  nonlinear  wave  conditions  on  certain  theoretical  and 
empirical  expressions.  Both  the  horizontal  and  vertical  components  of  velocity  are 
analysed.  Surprisingly  good  results  are  obtained  for  vertical  velocities,  even  for 
measurements  ne?’-  the  bottom  where  one  would  expect  difficulty  in  recording  the 
relatively  small  velocity  fluctuations. 

An  array  of  three  surface  piercing  capacitance  type  wave  gauges  was  placed 
in  the  middle  of  the  tank  to  resolve  incident  and  reflected  wave  height  spectra 
using  the  method  of  Goda  and  Suzuki  (1976).  This  method  is  itself  based  on 
linear  superposition  of  individual  wave  components.  However,  the  reflection  coef¬ 
ficients  were  on  the  order  of  5-15%  thereby  minimizing  the  importance  of  a  precise 
resolution  of  incident  and  reflected  spectra. 

Before  the  irregular  wave  time  series,  regular  waves  were  run  in  order  to  op¬ 
timize  the  design  of  the  wave  absorbing  beach.  It  was  discovered  that  the  reflec¬ 
tion  coefficient  varies  directly  with  the  ratio  wave  height/depth.  Surprisingly,  the 
amount  of  wave  reflection  was  not  correlated  with  peak  spectral  period. 

A  total  of  108  irregular  wave  time  series  were  run.  The  water  depth,  peak 
spectral  period,  and  significant  wave  height  were  varied  to  provide  a  range  of  wave 
conditions.  LDV  measurements  were  made  at  mid-depth  and  at  3  cm  above  the 
bottom.  Each  time  series  was  run  for  approximately  5  min.  Collection  of  velocity 
data  with  the  LDV  and  water  elevation  data  with  the  capacitance  gages  was  not 
synoptic.  Instead,  statistics  based  on  measured  velocity  data  were  compared  with 
theoretical  statistics  based  on  the  measured  water  elevations. 

One  objective  of  this  study  is  the  investigation  of  theoretical  expressions  for  the 
root -mean- square  horizontal  and  vertical  velocity  (Equations  4  and  5).  The  rms 
velocity  is  a  useful  and  important  parameter  in  studies  of  sediment  transport  and 
turbulence.  Results  indicate  that  predictions  within  10%  are  possible  given  the 
components  of  the  incident  wave  height  spectrum,  the  reflection  coefficients,  and 
the  water  depth.  Empirical  corrections  to  Equations  4  and  5  are  not  developed 
because  of  the  lack  of  clear  trends  between  the  measured  and  theoretical  values. 
One  would  expect  only  a  small  deviation  between  predicted  and  measured  urrrl, 
and  wrm_,  due  to  nonlinearities.  This  is  due  to  the  fact  that  while  the  nonlinear 
wave  profile  is  different,  the  value  of  rjrm>  is  only  slightly  different.  The  same  effect 
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would  extend  to  velocities. 


In  applying  Equations  4  and  5  this  investigation  discovered  that  there  is  no 
difference  between  using  a  random  phase  angle  (0n)  and  a  measured  phase  angle. 
A  measured  phase  angle  results  from  the  Fourier  analysis  of  the  wave  signal.  A 
random  phase  angle  is  a  random  value  between  0°  and  90°  applied  in  Equations  4 
and  5  to  predict  urTn,  and  wTTn,.  Also,  the  range  of  variability  of  the  rms  velocities 
with  varying  distance  from  the  toe  of  slope  (x)  is  less  than  0.05%  for  horizontal 
distances  more  than  1  m  from  the  reflecting  structure  (This  result  may  not  be 
applicable  to  structures  with  reflection  coefficients  greater  than  20%).  Therefore, 
it  is  concluded  that  the  value  of  x  in  Equations  4  and  5  should  be  simply  the 
distance  from  the  toe  of  slope. 

Also  studied  is  the  application  of  an  extended  form  of  the  Gaussian  probability 
distribution  function  to  wave  particle  velocities  (Equation  6).  That  is,  the  distri¬ 
bution  of  all  velocity  measurements,  not  just  peak  velocities.  This  Gram-Charlier 
distribution  is  in  the  form  of  a  series  using  Hermite  polynomials  and  makes  use  of 
higher  moments  of  the  distribution  like  skewness  and  kurtosis.  It  is  concluded  that 
the  theoretical  Gram-Charlier  distribution  is  less  effective  than  the  plain  Gaussian 
distribution  in  fitting  the  measured  data. 

Analysis  of  peak  velocities  reveals  that  a  modified  form  of  the  Rayleigh  prob¬ 
ability  distribution  termed  the  “Beta- Rayleigh”  is  an  improvement  over  the  plain 
Rayleigh  distribution  in  comparison  with  measured  data.  The  degree  of  improve¬ 
ment  varies  depending  on  velocity  direction  and  location.  The  results  are  summa¬ 
rized  in  Table  4.  The  Beta-Rayleigh  distribution  is  a  three  parameter  distribution 
{Umax,  UTma,  UTm q)  while  the  Rayleigh  distribution  is  only  one  parameter.  Another 
impiovement  is  that  the  Beta-Rayleigh  distribution  is  bounded  by  a  maximum  ve¬ 
locity  whereas  the  Rayleigh  distribution  “tail”  trails  off  into  infinity.  Also,  the 
Beta-Rayleigh  distribution  transforms  itself  into  the  plain  Rayleigh  distribution 
when  Umax  approaches  infinity  and  UTmq  —  2y/2 U?mt. 

Empirical  equations  are  developed  to  estimate  the  parameters  of  the  Beta- 
Rayleigh  distribution  in  terms  of  the  following  measured  parameters  of  the  incident 
wave  spectrum;  Tp,  Hmo,  and  water  depth,  h.  Good  results  are  achieved  with  the 
empirical  equations.  Correlation  coefficients  varied  between  0.78  and  0.97.  Fur¬ 
ther  research  is  recommended  to  broaden  the  data  base  from  which  the  empirical 
equations  are  developed.  Field  studies  would  be  particularly  helpful.  Additional 
efforts  are  recommended  to  investigate  the  change  with  elevation  (r)  of  the  equa¬ 
tions.  Eventually,  equations  could  be  developed  that  are  a  function  of  z,  and  are 
applicable  throughout  the  water  column. 
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Additional  recommendations  include  investigations  to  determine  the  variability 
of  velocity  statistics  with  depth,  z,  as  well  as  investigations  to  determine  the 
thickness  of  the  boundary  layer  under  oscillatory  flow.  Future  studies  also  could 
look  at  theoretical  wave  reflection  from  a  variable  slope  and  slopes  where  the 
reflection  coefficient  is  greater  than  20%. 
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APPENDIX  I:  SUPPLEMENTAL  GRAPHS 
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Figure  38:  uTTn,  Measured/Predicted  vs.  Spectral  Significant  Wave  Height 


urms(measured)/urms(predicted) 
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Figure  39:  uTm,  Measured/Predicted  vs.  Water  Depth 


s(measured)/urms(predicted) 


Measured/Predicted  vs. 


urms(measured)/urms(predicted) 
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Figure  42:  urm.  Measured/ Predicted  vs. 


wrms(measured)/wrms(predicted) 


Measured/ Predicted  vs.  Peak  Spectral  Period 


wrms(measured)/wrms(predicted) 
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Figure  45:  wrm,  Measured/Predicted  vs.  Spectral  Significant  Wave  Height 


wrms(measured)/wrms(predicted) 
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Figure  46:  wrm,  Measured/Predicted  vs.  Water  Depth 
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Figure  47:  wrm,  Measured/Predicted  vs.  Z/h 


ms(measured)/wrms(predicted) 


are  48:  wTrn,  Measured/Predicted 


wrms(measured)/wrms(predicted) 


Figure  49:  wTTna  Measured/Predicted  vs.  Ursell  No.  (without  outliers) 


wrms(measured)/wrms(predicted) 


Figure  50:  wTm,  Measured/Predicted  vs.  ZHmo/h 
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Figure  51:  %-Range  of  urm.  Predicted  as  a  function  of  X,  (Measured  Phase  Angle) 


e  53:  %-Range  of  wrmt  Predicted  as  a  function  of  X,  (Measured  Phase  Angle) 
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Figure  54:  No.  Waves  Counted,  Horizontal  Direction,  (a)  vs.  (b)  vs. 

Ursell  No. 
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Figure  55:  No.  Waves  Counted,  Vertical  Direction,  (a)  vs.  Hmo/h ;  (b)  vs.  Ursell 
No. 
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Figure  56:  %-Root-mean-square  Error,  Beta-Rayleigh  Distribution,  Positive  Hor¬ 
izontal  Direction,  (a)  Mid-depth;  (b)  Near  Bottom 
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Figure  57:  %- Root-mean-square  Error,  Beta-Rayleigh  Distribution,  Positive  Ver¬ 
tical  Direction,  (a)  Mid-depth;  (b)  Near  Bottom 
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Figure  58:  %-Root-mean-square  Error,  Beta- Rayleigh  Distribution,  Negative  Ver¬ 
tical  Direction,  (a)  Mid-depth;  (b)  Near  Bottom 
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figure  59:  lhe  Lftect  Of  Varying  Filter  Width,  Horizontal  Velocity  Peaks,  Time 
Series  No.  170.  (a)  Raw  Data;  (b)  Filtered  Data  [0.20-5.00]  s;  (c)  Filtered  Data 
[0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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. j 


Series  No.  19.  (a)  Raw  Data;  (b)  Filtered  Data  [0.20-5.00]  s;  (c)  Filtered  Data 
[0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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Figure  61:  The  Effect  of  Varying  Filter  Width,  Horizcntal  Velocity  Peaks,  Time 
Series  No.  70.  (a)  Raw  Data;  (b)  Filtered  Data  [0.20-5.00]  s;  (c)  Filtered  Data 
[0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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Figure  62:  The  Effect  of  Varying  Filter  Width,  Vertical  Velocity  Peaks,  Time 
Series  No.  170.  (a)  Raw  Data;  (b)  Filtered  Data  fO. 20-5. 00]  s;  (c)  Filtered  Data 
[0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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Figure  63:  The  Effect  of  Varying  Filter  Width,  Vertical  Velocity  Peaks,  Time 
Series  No.  19.  (a)  Raw  Data;  (b)  Filtered  Data  [0.20-5.00]  s;  (c)  Filtered  Data 
[0.20-10.00]  s;  (d)  Filtered  Data  [0.10-5.00]  s 
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Figure  66:  Maximum  Velocity  Measured/Predicted  vs.  Relative  Depth,  (a)  Hori¬ 
zontal  Direction;  (b)  Vertical  Direction 
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Figure  69:  Parameterization  of  Beta- Rayleigh  Distribution,  Urm »•  (a)  VS.  Hmo/h 
(b)  vs.  Urseli  No. 
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Figure  70:  Parameterization  of  Beta- Rayleigh  Distribution,  f7rmq.  (a)  vs.  Hmo/h 
(b)  vs.  Ursell  No. 
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Figure  71:  Parameterization  of  Beta- Rayleigh  Distribution,  WTmt.  (a)  vs.  Hmo/h ; 
(b)  vs.  Ursell  No. 
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Figure  72:  Parameterization  of  Beta- Rayleigh  Distribution,  Wrmq.  (a)  vs.  Hmo/h 
(b)  vs.  Ursell  No. 
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Figure  73:  Time  Series  No.  9.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Figure  74:  Time  Series  No.  153.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0  20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  75:  Time  Series  No.  18.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Fntered  Data 
[0  °0-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P(U)  Prob.  Density  -  P(W) 
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Figure  76:  Time  Series  No.  19.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  77:  Time  Series  No.  25.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P<U)  Ptob.  Density  -  P(W) 
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Figure  78:  Time  Series  No.  28.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P(U)  __  Prob.  Density  -  P(W) 
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Figure  79:  Time  Series  No.  29.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  80:  Time  Series  No.  35.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P(U)  Prob.  Density  -  P(W) 
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Figure  81:  Time  Series  No.  38.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0. 20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P(W) 
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Figure  82:  Time  Series  No.  39.  (a)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  83:  Time  Series  No.  45.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  84:  Time  Series  No.  48.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Proto.  Oanaily  -  P(U)  _  Prob.  Density  -  P(W) 
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Figure  85:  Time  Series  No.  49.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Oanaity  -  P(U)  _  „  „  PfOb.  Dbnlty  -  P(W) 
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Figure  86:  Time  Series  No.  55.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


_  Proto.  0«n«ty  -  P(U)  Proto.  Density  -  P(W) 

0.26  | - - - 1  0.12  i - 
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Figure  87:  Time  Series  No.  58.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0. 20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Proto.  Dontity  -  P(U)  _  .  Proto.  Donoity  -  P(W) 

, - - — -  - 1  0.12  i - 
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Figure  89:  Time  Series  No.  65.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0  20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Dcntity  -  PON) 
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Figure  90:  Time  Series  No.  68.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P<U)  Prob.  Density  -  P(W) 
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Figure  91:  Time  Series  No.  69.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Density  -  P<U)  Pfob.  Density  -  P(W) 
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Figure  92:  Time  Series  No.  75.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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Figure  93:  Time  Series  No.  78.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Den»ity  -  P(U)  _  _  Prob.  Den#ity  -  P(W) 
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Figure  94:  Time  Series  No.  79.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


Prob.  Osnsity  -  P(U)  _  .  PfQb.  Den«ity  -  P(W) 
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Figure  95:  Time  Series  No.  85.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 


PfQfa.  Oantily  -  P(U) _ _ _  0  , 2  Pf°b-  0<>n,ity  ~  P(W) 
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Figure  96:  Time  Series  No.  88.  (a)  Filtered  Data  [0.20-5.00]  s;  (b)  Filtered  Data 
[0.20-5.00]  s;  (c)  Raw  Data;  (d)  Raw  Data 
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amplitude  of  incident  wave; 
amplitude  of  nth  linear  component; 
frequency; 

gravitational  acceleration; 
water  depth; 

monochromatic  wave  height; 
breaking  wave  height; 

zeroth-moment  wave  height  equal  to  four  times  the  standard 

deviation  of  sea  surface  elevations; 

root-mean-squared  wave  height  of  irregular  waves; 

nth  order  Hermetian  polynomial; 

wavenumber  of  monochromatic  wave  =  2tt /  i ; 

wavenumber  of  the  nth  linear  component; 

reflection  coefficient  for  a  monochromatic  wave; 

reflection  coefficient  of  the  nth  lineai  component; 

monochromatic  wave  length; 

deep  water  monochromatic  wave  length; 

time; 

wave  period; 
mean  wave  period; 

wave  period  associated  with  the  peak  spectral  frequency; 

horizontal  component  of  fluid  velocity; 

root -mean- square  of  horizontal  velocity  time  series; 

maximum  horizontal  velocity  over  a  wave  cycle; 

mean  of  the  horizontal  velocity  peaks; 

root-mear.  of  the  horizontal  velocity  peaks; 

root-mean-quad  of  the  horizontal  velocity  peaks; 

vertical  component  of  fluid  velocity; 

root-mean-square  of  vertical  velocity  time  series; 

maximum  vertical  velocity  over  a  wave  cycle; 

mean  of  the  vertical  velocity  peaks; 

root-mean-square  of  the  vertical  velocity  peaks; 

root-mean-quad  of  the  vertical  velocity  peaks; 

horizontal  dimension; 

vertical  dimension; 

parameter  of  the  Beta-Rayleigh  distribution; 
parameter  of  the  Beta-Rayleigh  distribution; 

Gamma  function; 
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sea  surface  elevation; 

root-mean-square  of  sea  surface  elevations; 
phase  shift  due  to  reflection; 

phase  shift  of  the  nth  linear  component  due  to  reflection; 

skewness; 

kurtosis; 

higher  order  moment  of  the  nondimensional  probability 
distribution  function; 
mathematical  Pi  =  3.14159. . . ; 
angular  wave  frequency; 

angular  wave  frequency  of  the  nth  linear  component. 
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